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Director’s  Foreward 

Severzd  months  ago  I  asked  a  group  of  6RL  scientists  and  engineers  to  put  together 
a  “White  Paper”  on  the  knowledge  gaps  that  most  seriously  deter  our  ability  for 
analytical  and  computationad  prediction  of  phenomena  in  armor,  penetration,  and 
warhead  mechanics.  My  purpose  was  to  gjun  focus  for  our  internal  efforts  auad  to 
express  our  needs  to  others  who  might  join  our  quest.  I  am  optimistic  that  others 
will  do  so  since  our  needs  are  not  unique.  I  am  sure  that  many  will  recognize 
that  we  w’ill  accomplish  more  by  cooperative,  focused  effort  then  we  could  do 
separately. 

This  document  is  the  result  of  my  request.  It  overwhelmingly  exceeds  my 
expectations  which  were  for  a  summary  list  with  a  brief  rationalization  for  each. 
It  provides  focus  for  specific,  priority  work  and  also  gives  a  broad  perspective  as 
context  for  these  needs. 

My  hope  is  that  this  will  be  a  living  document  that  is  nurtured  by  the  BRL  and 
by  the  community  at  large  confronted  with  problems  involving  material  response 
at  extremes  of  load  and  deformation  and  their  rates.  I  am  certain  it  will  prove 
invaluable  if  it  helps  to  maintain  an  evolving  focus  on  our  collective  needs.  I  invite 
everyone  involved  in  the  challenging  problems  represented  here  to  join  with  us  in 
defining  and  maintaining  this  focus. 
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Preface  (T,  Wright) 

Historical  Remarks  Just  over  a  decade  ago  the  National  Materials  Advisory 
Board  (NMAB)  published  a  report  entitled  “Matericils  Response  to  Ultra-High 
Loading  Rates”  (1).  This  report  surveyed  the  state  of  material  modeling  then  in 
existence,  is  appropriate  for  analysis  and  design  of  conventional  ordn<uice  devices. 
By  way  of  background,  NMAB-356  came  into  existence  because  it  was  widely 
recognized  by  the  mid  1970’s  that  inadequate  description  of  severe  deformation 
and  failure  in  solids  constituted  a  major  limitation  on  the  effectiveness  of  large 
scale  computations  in  ordnance  applications.  Three  years  prior  to  publication 
of  the  report,  at  the  urging  of  the  B^dlistic  Research  Laboratory  (BRL)  and  the 
Army  Research  Office  (ARO),  the  office  of  the  Director  of  Defense  Research  and 
Engineering  (DDR&E,  DOD)  had  requested  NMAB  to  assess  the  state  of  the  art 
as  it  existed  at  that  time,  and  to  recommend  future  directions  for  research  and 
development.  Consequently  an  expert  committee,  drawn  from  academia  and  from 
both  private  and  government  reseaich  laboratories,  was  formed,  and  cifter  a  series 
of  meetings,  NMAB-356  was  born. 

For  various  reasons  the  committee  confined  itself  to  a  discussion  of  metals  only. 
At  the  risk  of  gross  oversimplification,  some  of  the  major  points  that  emerged  may 
be  summarized  as  follows.  In  decreasing  order  of  importance,  limitations  in  the 
state  of  material  modeling  were  found  to  be:  __  _ 

1.  severe  in  dynamic  material  failure; 

2.  moderate  in  plastic  deformation; 

3.  trivial  for  most  purposes  in  equations  of  state. 

Recommended  research  emphasized:  nucleation,  growth,  and  coalescence  of  voids  • 
in  ductile  materials  and  of  cracks  in  brittle  materials;  edl  aspects  of  the  formation 
and  propagation  of  adiabatic  shear  bands;  and  development  of  the  underlying 
micromechanics,  thermodynamics,  and  kinematics  of  plastic  deformation  within 
the  context,  of  a  consistent,  finite  deformation  theory.  The  report  also  suggested 
the  development  of  new  test  methods  for  high  rate  deformation  and  fmlure.  The 
interested  reader  is  referred  to  the  originad  report  for  many  fascinating  details  and 
partiad  analyses,  which  in  toto  give  a  vivid  snapshot  qf  the  state  of  the  art  in  dy¬ 
namic  materiad  modeling  circa  1980.  NMAB-356  served  ais  a  cat^yst  and  rationale, 
for  much  new  work  that  was  initiated  in  the  early  years  of  the  decade  following 
its  publication,  particularly  in  the  area  of  adiabatic  shear  bands.  Predictably, 
research  in  the  academic  and  defense  communities  did  not  develop  precisely  along 
the  lines  called  tor,  but  on  the  whole  NMAB-356  was  a  remarkably  prescient 
document. 

Purpose  T:)day,  twelve  years  later,  as  the  tremendous  advtinces  in  comput¬ 
ing  power  made  during  the  1980’s  continue  to  accelerate,  it  is  feasible  to  devote 
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an  increasing  fraction  of  the  computing  cycle  to  the  evaluation  of  more  sophis¬ 
ticated  constitutive  laws.  Therefore,  it  seems  appropriate  to  reassess  otir  needs 
in  the  area  of  materia  modeling  for  terminal  b2Lllistic  computation  and  analysis. 
Toward  this  end,  this  report  presents  an  overview  of  the  current  state  of  material 
modeling  in  terminal  ballistics  codes  and,  perhaps  more  importantly,  describes 
th<^  deficiencies  in  the  state  of  material  modeling  which  are  perceived  to  impose 
the  most  severe  limitations  on  the  capability  amd  reliability  of  terminal  ballistic 
simulation.  Conversely,  this  report  is  not  intended  to  be  a  comprehensive  treatise 
on  either  terming  ballistics  or  the  dynamic  response  of  solids  in  general.  The 
point  of  view  expressed  is  that  of  a  collection  of  professionals  within  the  Terminal 
Ballistics  Division  of  the  Ballistic  Research  Laboratory,  and  therefore  is  by  no 
means  comprehensive. 

User’s  Guide  The  report  opens  with  a  discussion  of  certain  aspects  of  ballis¬ 
tic  penetration  from  a  phenomenological  point  of  view.  Extensive  ballistic  testing 
has  demonstrated  that  adiabatic  shear  banding  is  a  major  mechanism  for  erosion 
''bfrnietals,  and  that  it  can  often  compete  with  bulk  plastic  flow  in  moving  material 
away  from  the  path  of  penetration.  Similarly  testing  has  demonstrated  that  tran¬ 
sient  effects  from  wave  propagation  must  always  be  considered  in  brittle  materials 
such  as  ceramics  and  cermets.  The  next  section  describes  some  of  the  major  codes, 
as  well  as  the  principal  material  models,  that  are  in  use  today  for  simulation  of 
terminal  effects.  The  next  four  sections  deal  with  specific  research  areas  that  seem 
to  have  particular  relevance  today.  Section  3  covers  in-house  and  related  efforts 
in  finite  plasticity,  section  4  describes  recent  results  in  modeling  adiabatic  shear 
bands,  section  5  is  coiKemed  with  brittle  failure  and  granular  flow,  and  section 
6  points  out  some  problems  with  equations  of  state  as  they  are  actually  used  in 
todays  impact  codes.  The  report  closes  with  a  summary  of  important  points  that 
have  arisen  in  the  body  of  the  report  and  a, statement  of  needs  that  address  those 
same  points. 

Closure  and  Call  for  Comments  Clearly  the  major  concerns  of  NMAfe- 
356  are  still  of  interest  today.  Failure  is  still  the  least  understood  aspect  of  material 
response,  and  therefore  is  receiving  the  most  attention.  However,  failure  in  metals 
preceded  by  finite  plastic  deformation,  which  determines  the  history  and  sets  the 
environment  within  which  damage  and  failure  evolve.  Analogously,  transient  wave 
propagation  creates  damage  in  brittle  materials  which  sets  the  environment  from 
which  subsequent  comminution  evolves.  Increased  use  of  cerzunics  hais  brought 
recognition  that  post  failure  response  is  also  important,  hence  the  interest  in 
granular  flows.  Finally,  it  seeriis  that  even  as  mature  a  subject  as  modeling  of 
equations  of  state  requires  some  care  and  attention. 

Interested  readers  are  invited  to  communicate  their  views  on  any  and  all  as¬ 
pects  of  this  report  to:  ' 
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Director 

BaUistic  Research  Laboratory 
Attn:  SLCBR-D 

Aberdeen  Proving  Ground,  Maryland  21005-5066 
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1  Phenomenology  of  Ballistic  Penetration 

\ 

14  Fundamentals  of  the  Penetration  Process  (L.  Magness) 

The  penetration  of  typiccil,  ductile,  metallic  armor  by  a  high-velocity  kinetic  en- 
ergy  penetrator  (a  rod-like  projectile)  or  by  the  still  higher  velocity  jet  of  metal 
produced  by  a  shaped  charge  wao'head  can,  to  the  first  order,  be  viewed  as  a  hy¬ 
drodynamic  process*,  as  depicted  in  Figure  1.  After  a  transient  phase  upon  ini  tied 
impact,  the  penetrator-target  interaction  quickly  settles  down  to  a  nearly  steady 
process  in  which  the  head  of  the  impacting  projectile  or  jet  is  eroded  away  at  the 
penetrator-target  interface,  while  the  interface  moves  forward  displacing  a  cavity 
in  the  armor.  The  kinetic  energy  of  the  impacting  projectile  or  jet  is  expended 
prim2iriiy  as  plastic  work,  performed  in  opening  a  cavity  in  the  armor  material 
2md  eroding  and  inverting  the  penetrator  itself.  Very  large  peak  stresses  and 
plastic  strain  rates  occur  during  the  initial  impact  phase,  but  these  levels  quickly 
diminish  to  lower  values  during  the  nauch  longer  steady  phase  of  the  penetration 
process  [2].  However,  due  to  the  deceleration  of  the  incoming  penetrator  material ' 
at  the  penetrator-target  interface  and  to  the  corresponding  acceleration  of  the  tar¬ 
get  material  out  of  the  path  of  the  moving  interface,  the  inertial  forces  developed 
during  the  steady-state  phase  will  still  grci  dy  exceed  the  yield  strengths  of  both 
the  target  and  penetrator  materials. 

First  order  analytical  penetration  models  treat  the  entire  penetration  process 
as  simply  a  fluid-on- fluid  interaction,  Birkhoff  [3j.  An  expression  for  the  final 
depth  of  penetration,  P,  derived  from  the  Bernoulli  equation,  results: 


where  pp  and  p,  are  the  mass  densities  of  the  penetrator  and  target  materials 
respectively,  and  L  is  the  length  of  the  penetrator.  This  often  gives  quite  reason¬ 
able  penetration  estimates  for  the  performance  of  shaped  charge  jets,  when  L  is 
corrected  for  stretching  of  the  jet.  With  the  inertial  forces  generated  by  the  high 
velocities  of  jets  dominating  the  penetration  process,  the  penetration  of  ductile 
metallic  targets  thus  becomes  largely  a  function  of  jet  density  and  length  alone. 
This  simple  expression  explciins  the  use  of  high-ductility  copper  in  shaped  charge 
liners,  to  obtain  stretching  jets  of  greater  lengths,  and  the  use  of  high-density 
materials  such  as  uranium  and  tantalum  for  shaped  charge  liners,  and  of  uranium 
or  tungsten  for  long-rod  penetrators. 

The  fluid  dynamic  model  is  less  applicable  for  impacts  at  somewhat  lower  veloc¬ 
ities,  such  as  the  ordnance  velocities  of  modern  fin-stabilized,  long  rod  penetrator 
munitions,  or  when  the  mechanical  responses  of  the  penetrator  or  target  during 

*It  should  be  emphasized  tha;  penetration  into  a  ceramic  target 'produces  a  very  different 
response  than  that  discussed  here;  see  Section  1.2 


the  penetration  process  axe  other  than  ductile  erosion  and  cavity  expamsion,  re¬ 
spectively.  Additional  terms  representing  the  yield  strengths  of  the  penetrator  and 
target  materials  are  often  included  in  modified  versions  of  the  Bernoulli  equation  to 
improve  the  representation  of  penetration  events  by  these  analytical  models.  The 
eroding-rod  models  (e.g.  Tate  [4,  5],  Alekseevski  [6|,  Frank  and  Zook  [7])  which  in¬ 
clude  these  strength  terms  reflect  the  experimentally  observed,  velocity-dependent 
penetration  of  long-rod  projectiles  (Figure  2)  and,  when  properly  calibrated,  can 
give  reasonably  accurate  predictions  of  their  overall  penetration  perfonAance. 

In  the  case  of  brittle  materials  (ceramics),  the  assumptions  of  the  steady-state, 
hydrodynamic  model  must  be  questioned  not  only  from  the  viewpoint  of  time- 
dependent  strength  degradation  but  also  becatise  experiments  indicate  a  major 
contribution  from  transient  wave  interaction  effects.  This  point  is  considered  in 
more  detail  in  the  next  section. 

In  reality,  neither  the  penetrator  nor  the  target  materials  behave  as  fluids  or 
modified  fluids.  To  better  understand  the  penetration  process,  finite  element  or 
finite  difference  numerical  simulations  are  employed.  Two-  and  three-dimensional 
hydrocodes,  both  Eulerian  and  Lagrangian,  attempt  to  model  the  actual  deforma¬ 
tion,  flow,  amd  fracture  behavior  of  both  penetr  ator  and  target  materials.  Through 
constitutive  equations  for  each  material,  their  mechanical  responses  under  the 
high  stress  and  high  strain  rate  loading  cor  lit’  .  s  are  calculated  and  followed 
through  each  time  increment  of  the  peuetrat’  .'r  j^rocess.  The  ability  of  the  nu¬ 
merical  models  to  simulate  the  penetration  £• .  '.it  thus  depends  on  the  accuracy 
and  completeness  of  the  materials  models  anc  o;  the  computational  algorithms 
employed  in  the  hydrocode  itself.  In  the  ne  ^t  section,  several  examples  of  some 
of  the  coiiiplexities  of  penetrator-target  iiitfiartions  which  have  been  observed 
experimentally  are  given.  Aspects  of  the  mt  chaiucal  responses  of  penetrator  and 
target  materials  that  have  been  shown  to  be  ‘mportant  to  the  understanding  and 
prediction  of  ballistic  test  results  are  illustr -.'.."d  by  two  phenomenological  exam¬ 
ples:  (1)  the  penetration  of  hard  ceramics,  and  (2)  penetrations  involving  ductile 
materiaJs,  both  armor  and  different  high  den. ^I'y  penetrator  materials.  In  the  sub¬ 
sequent  sections,  the  ability  of  the  computational  hydrocodes  to  model  these  and 
other  aspects  of  the  penetration  process  will  be  reviewed  and  critiqued. 

1.2  Phenomenological  Observations  of  Penetration  Into 
Brittle  Materials  (G.  Hauver) 

Ceramics  are  commonly  assumed  to  have  an  inherent  resistance  to  penetration 
by  a  long-rod  projectile.  However,  test  results  often  provide  evidence  that  the 
ballistic  performance  of  a  ceramic-laminate  target  varies  with  time  during  the 
course  of  penetration.  For  example,  increa-sing  the  proportion  of  ceramic  in  a 
target  commonly  fails  to  produce  a  corresponding  improvement  in  its  resistance 


to  penetration.  In  recent  studies  at  the  Ballistic  Research  Laboratory  (BRL), 
all  available  data  from  quarter-scale  ballistic  tests  were  analyzed  with  the  aid 
of  a  penetration  model  [4,  5]  using  the  target  resistance  term,  i2,  as  a  measure 
of  ceramic  performance.  Such  steady-state  models  are  not  intended  to  describe 
penetration  when  the  resistance  term  is  time  dependent.  However,  the  steady- 
state  condition  was  satisfied  by  using  a  point-by-point  analysis  and  assuming 
an  average  value  of  resistance  from  the  beginning  of  penetration  to  a  specific 
penetration  depth  in  the  ceramic  of  each  test.  This  approach  greatly  suppressed 
the  actual  change  in  resistance  during  penetration,  but  it  preserved  performance 
trends  which  could  be  related  to  sources  of  damage  to  the  ceramic  that  cause  its 
R  to  ch2uige.  Supporting  tests  were  then  conducted  to  identify  features  of  target 
design  which  influenced  the  damage.  - 

The  analysis  of  quarter-scale  bcdlistic  tests  indicated  similar  behavior  for  dif¬ 
ferent  ceramics,  and  representative  performance  trends  are  shown  in  Figure  3.  In 
this  figure,  average  resistance  to  penetration,  R,  is  plotted  as  a  function  of  time, 
T.  With  no  cover  plate  at  the  front  of  the  ceramic,  performance  typically  followed 
Curve  A.  The  performance  was  initially  low,  but  it  increased  with  time,  passed 
thrpugh  a  maximum,  and  decreased  at  later  times  when  different  sources  of  dcim- 
age  became  dominant.  However,  the  decrease  in  performance  at  late  times  did  not 
always  occur.  In  one  case  [8]  the  performance  continued  to  increase  along  Curve 
AA.  This  trend  occurred  when  the  ceramic,  without  a  cover  plate,  was  heavily 
confined  at  the  back  and  side.  With  a  steel  cover  plate,  performance  typically 
followed  Curve  B.  The  presence  of  a  cover  plate  suppressed  early  damage,  but 
the  curve  plunged  rapidly  as  other  damage  sources  influenced  the  perfotmance. 
However,  with  efforts  to  suppress  damage  in  tests  at  the  BRL,  performance  of 
the  same  ceramic  could  follow  a  higher  curve,  BB.  At  present,  maximum  perfor¬ 
mance  (highest  position  of  Curve  BB)  has  not  been  determined.  It  is  known  that 
if  damage  is  not  suppressed  and  performance  follows  Curves  A  or  B,  the  late-time 
performance  of  most  ceramics  lies  within  a  narrow  range  of  resistance  values.  It 
is  this  late-time  performance  which  is  often  used  to  rank  the  ballistic  capability 
of  ceramics. 

Tau'get  performance  depends  on  damage  to  the  ceramic  component.  Although 
this  damage  is  influenced  by  characteristics  of  the  ceramic,  it  is  strongly  depen¬ 
dent  on  response  of  the  total  target  system.  Ceramics,  in, general,  are  strong  in 
compression  but  weak  in  tension.  Although  high-quality  cersimics  may  respond 
elastically  to  compressive  stresses  above  15  GPa  in  uniaxial  strain,  they  may  fail 
at  tensile  stresses  below  0.5  GPa  [9,  10].  The  design  of  targets  for  ballistic  evalu¬ 
ations  of  ceramics  or  of  actual  armor  packages  should  minimize  conditions  which 
produce  damage  and  provide  an  environment  in  which  the  compressive  strength 
of  the  ceramic  may  be  exploited  for  ballistic  protection.  The  damage  produced 
by  a  long-rod  penetrator  is  found  to  depend  on  many  aspects  of  target  design  and 
response.  Current  information  about  target  behavior  and  sources  of  damage  is 


briefly  siunmarized  in  tbe  following  paragraphs. 

Impact  Damage  Tests  with  and  without  a  steel  plate  in  front  of  the  ceramic 
re^^al  a  large  difference  in  the  r^istance  to  penetration,  especially  at  early  tipies 
after  imp^.  A  cover  plate  is  believed  to  offer  two  benefits.  First,  it  prevents 
the  cer^nic  from  experiencing  the  peak  stress  which  may  be  40  to  50  GPa  during 
rect  impact  by  a  high-density  penetrator  traveling  at  current  ordnance  velocities. 
Second,  a  steel  ojver  plate  provides  confinement  which  retards  the  displacement 
of  cerwiic  materiad,  as  a  block,  from  the  front  of  the  target;  it  also  suppresses  an 
almost  explosive  failure  which  would  occur  if  the  highest  stresses  were  relieved  at 
a  free  boundary  of  the  ceramic. 

Confinement  Structural  support  must  also  be  provided  at  the  back  and 
side  of  a  c«!ramic.*  When  a  long-rod  penetrator  arrives  at  the  ceramic,  radial 
cracks  quickly  develop,  producing  wedge-shaped  blocks  w'hich  load  the  side  con¬ 
finement.  Weakly  constructed  confinement  may  fail  catastrophically.  With  more 
substantial  confinement,  the  surrounding  ductile  material  is  engraved  by  each 
wedge.  Collectively,  the  radial  motion  of  the  wedge-shaped  blocks  increases  the 
total  cross-sectional  area  of  the  target  occupied  by  the  ceramic.  Even  with  what 
may  be  referred  to  as  heavy  confinement,  the  ceramic  cross-sectional  area  may 
increase  by  several  times  the  cross-sectional  area  of  the  penetrator.  Consequently, 
the  interaction  of  the  penetrator’s  lateral  surface  with  the  penetration  “tunner 
produced  in  the  target  may  be  much  weaker  than  it  would  be  if  radial  displace¬ 
ment  of  the  wedges  did  not  occur.  As  an  example.  Figure  4  shows  the  increase 
in  area  determined  during  penetration  into  a  target  used  for  ballistic  evaluations. 
Penetration  through  the  steel  cover  was  complete  at  time  Tl  and  penetration 
reawdied  the  steel  backing  at  time  T2.  By  time  T2  the  total  cross-sectional  area 
occupied  by  the  ceramic  had  increa.sed  by  nearly  11  times  the  cross-sectional  area 
of  the  penetrator.  Extrapolation  to  the  abscissa  suggested  an  early  development 
of  ceramic  wedges.  Deviation  of  the  actual  curve  indicated  confinement  motion 
resulting  from  arrival  of  the  impact  shock.  Internal  stresses  also  act  during  pen¬ 
etration  to  displace  confinement  at  the  front  and  back  of  the  ceramic.  The  cover 
plate  tends  to  bulge  outward  emd  interfaces  of  a  layered  target  tend  to  separate, 
even  if  the  deformation  of  clamping  structures  remains  within  the  ela.stic  range. 
If  interfaces  extend  directly  to  the  side  boundary  of  the  target,  erosion  products 
from  the  penetrator  and  small  particles  of  ceramic  are  vented  to  the  exterior. 

Wave  Propagation  and  Interaction  Most  efforts  to  confine  the  ceramic 
component  tend  to  consider  structural  support  while  neglecting  wave  propaga¬ 
tion  in  the  finite  target.  A  number  of  factors  influence  wave  propagation  and 
interaction.  ' 

1.  A  low-impedance  bonding  layer  is  usually  present  between  ceramic  tile-^  and 
between  the  ceramic  and  its  confinement.  This  layer  is  a  source  of  tensile  le 
lease  and  associated  damage  to  the  ceramic.  Thickness  of  the  interface  layer 


may  not  be  well  controlled,  and  air  may  be  entrapped,  further  contributing 
to  variable  dfimage  and  performance. 

2.  Even  if  the  interface  and  confinement  could  be  perfectly  matched  to  the 
acoustic  impedance  of  the  ceramic,  confinement  thicknesses  <ire  usucdly  fi¬ 
nite  and  tensile  reflections  soon  return  from  the  free  boundary.  Wave  re¬ 
flections  within  the  ceramic  then  depend  on  the  gwmetry  of  the  reflecting 
free  boundary.  Tests  have  shown  that  the  damage  and  loss  of  performance 
caused  by  the  interaction  of  high-amplitude  reflected  waves  Cein  be  reduced 
by  installing  wave  (momentum)  traps  at  the  boundary  of  the  confinement. 
Figure  5  shows  orthogonal  post-test  radiographs  which  depict  plcinar  fciil- 
ure  produced  by  interacting  tensile  waves  reflected  from  the  boundary  of 
a  square  target.  Figure  6  shows  orthogonal  post-test  radiographs  which 
confirm  that  the  penetration  path  regains  symmetry  when  wave  traps  are 
added  at  the  free  boundary  of  a  square  target.  Figure  7  shows  fully  de¬ 
veloped  orthogonal  failure  planes  revealed  when  the  steel  cover  plate  was 
removed  from  a  square  target  without  wave  traps.  This  figure  also  shows 

•  the  wedge-shaped  blocks  through  which  force  at  the  penetration  front  is 
transmitted  to  the  side  confinement. 

3.  Test  results  suggest  that  low-amplitude  reflected  waves  may  contribute  to 
the  long-term  cumulative  damage  in  a  ceramic  target,  although  this  is  not 
well  established.  Wave  traps  lose  effectiveness  as  the  wave  amplitude  de¬ 
creases,  so  the  reflection  and  interaction  of  weak  stress  waves  is  not  easily 
controlled. 

Resistance  to  Penetration  Many  aspects  of  design  and  behavior  are  be¬ 
lieved  to  influence  the  resistance  of  a  ceramic  target  to  penetration  by  a  long  rod. 
Cover  plates  are  necessary  to  prevent  severe  impact  damage  to  the  ceramic,  but  rel¬ 
atively  slow  penetration  through  the  cover  plate  allows  time  for  the  faster  impact 
shock  to  damage  the  ceramic  before  the  penetrator  arrives.  Damage  by  intersect¬ 
ing  tensile  waves  from  a  square  target  boundwy,  and  the  resulting  degradation 
of  performance,  were  noted  earlier.  Converging  tensile  reflections  from  a  circular 
target  boundary  may  cause  a  greater  decrease  in  the  resistance  to  penetration. 
When  the  penetrator  arrives  at  ceramic  under  a  cover  plate,  it  conunonly  finds  a 
material  which  still  is  highly  resistant  to  penetration.  The  erosion  rate  increa.ses 
and  erosion  products  commonly  flow  into  the  interface  betw^n  the  cover  plate 
and  the  ceramic.  The  cover  plate  bulges  to  accommodate  the  erosion  products, 
while  the  ceramic  experiences  a  compressive  stress  which  may  suppress  damage 
and  increase  the  resistance  to  penetration.  Frequently,  t^lrget  construction  allows 
separation  of  the  entire  interface,  allowing  a  lateral  flow  of  erosion  products  and 
failed  ceramic  directly  to  the  outside  of  the  target,  relieving  internal  stresses.  This 
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behavior  may  promote  damage  and  reduce  the  resistance  of  the  ceramic  to  pene¬ 
tration.  As  penetration  proc^ds  into  the  ceramic,  major  failure  planes  develop, 
^^ng  blocks  of  cermnic  through  which  forces  at  the  penetration  front  are  coupled 
to  the  confinement.  At  the  side  of  the  penetration  path,  radial  cracks  produce 
pie  blocks  which  we  forced  against  the  side  confinement.  As  noted  before,  the 
increase  in  area  at  the  boundary  can  be  many  time  the  cross-section  J  area  of  the 
penetrator.  This  may  influence  granular  flow  in  the  penetration  path,  but  it  also 
may  provide  an  alternative  means  for  displacing  ceramic  from  in  front  of  the  pen¬ 
etrator.  Major  failure  planes,  intersecting  the  penetration  axis,  further  improve 
the  coupling  to  confinement  at  the  front  cind  back  of  the  ceramic.  This  must  con¬ 
tribute  to  deformation,  separation  of  interfaces,  additional  loss  of  ceramic  from 
the  target,  and  a  further  decrease  in  resistapce  to  penetration.  Sequential  failure 
of  the  confinement  may  make  a  greater  contribution  to  cumulative  damage  than 
long-term  wave  interactions.  Experimental  studies  may  soon  provide  an  answer. 

Summary  It  is  clear  that  the  mechanics  of  penetration  for  an  armor  contain¬ 
ing  ceramic  elements  is  very  different  frorri  the  purely  metallic  case.  The  concept 
of  intrinsic  resistance  is,  for  ceramic  elements,  almost  completely  spurious.  It 
has  become  evident  that,  due  to  the  mechanisms  just  discussed,  the  penetration 
resistance  of  a  ceramic  target  depends  fundamentcJly  on  the  response  of  the  total 
target  system.  Current  studies  at  the  BRL  show  that  modified  tcurget  configu¬ 
rations  can  produce  a  significant  increase  in  the  resistance  to  penetration,  but 
the  full  potential  for  improvement  has  not  yet  been  determined.  Current  ballistic 
twts  clearly  do  not  provide  an  environment  in  which  the  constitutive  response 
of  an  individual  ceramic  material  has  a  major  influence.  Indeed,  a  crucial  chal¬ 
lenge  in  the  analysis  of  ceramic  armors  is  to  understand  the  basic  penetration 
mechanics  well  enough  so  that  armor  packages  may  be  designed  which  do  exploit 
the  behaviors  of  different  ceramics.  It  is  further  apparent  that  to  obtain  reliable 
predictions  of  ballistic  performamce,  computational  siinulations  must  employ  more 
accurate  modeling  of  dynamic  material  behavior  than  is  now  the  case  and  must 
treat  the  response  of  the  total  target  system.  Substemtial  progress  will  proceed 
from  accurate  simulation  of  carefully  controlled  experiments. 


1.3  Phenomenological  Observations  of  Penetration  Involv¬ 
ing  Ductile  Materials  (L.  Magness,  M.  Raftenberg,  J.  Walter) 

While  the  hydrodynamic  picture  of  ballistic  penetration  (Figure  1)  is  appealing 
both  for  its  simplicity  and  successes,  it  ignores  an  essential  aspect  of  the  solid 
state:  internal  structure.  In  mechanical  terms  the  internal  structure  (microstruc¬ 
ture)  of  solids  is  expressed  as  material  anisotropy  and  inhomogeneity.'*  It  has 

*Bulk  specimens  of  metals  are  usually  aggregates  composed  of  many  crystals  which  are  always 
anisotropic  and  may  be  individually  inhomogeneous  doe  to  second  phases,  inclusions,  etc.  The 
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deep  penetration  tunnels  in  high-strength  steel  monoliths  reveals  shear  bands  that 
snwate  away  from  the  tunnel  in  a  spiral  pattern.  The  visible  bands  axe  the  trace 
in  the!  section  plame  of  two  families  of  (apparently)  counter-winding,  helical,  shear 
surf,^^  which  intersect  to  produce  pyramidal  fragments.  Often,  a  tube  composed 
of  these  fragments  (welded  together  because  of  the  high  temperatures  which  occur 
during  formation  of  the  shear  bands)  surrounds  the  penetrator  tunnel,  which  has  a 
roughened  wall  due  to  ejection  of  some  of  the  fragments.  Neither  feature  appears 
in  tMgets  of  lower  strength  steel. 

As  suggested  above,  distributed  shear  banding  tends  to  occur  at  higher  impact 
speeds  and  in  regions  (e.g.  near  the  target-penetrator  interface)  where  the  macro¬ 
scopic  rate  of  deformation  is  large  enough  to  initiate  band  growth  from  preexist¬ 
ing  or  deformation-induced  material  inhomogeneities.  Although  not  as  obvious  as 
plug-boundary  formation,  distributed  shear  band  formation  is  no  less  important 
foi  the  penetration  process.  In  this  mode  large  numbers  of  bcinds  nucleate  in  a 
more  or  less  dense  distribution  within  a  region  of  material.  As  the  bands  grow  the 
material’s  macroscopic  strength  is  reduced  anisotropically  because  the  bands  tend 
to  be  oriented  in  the  local  directions  of  maximum  Shear  stress  pr  shear  strain  rate. 
Here  again  the  target  presumably  offers  less  resistance  to  the  incoming  penetrator 
than  it  would  if  shear  bands  had  not  formed.  Ultimately,  the  distributed  shear 
bands  may  coalesce  to  form  fragments  of  the  target  material.  Upon  perforation, 
these  fragments  may  comprise  most  of  the  behind-armor  debris  which  is  one  of 
the  principal  lethal  mechanisms  of  the  penetrator.  Experimental  results  of  Curran 
et  al.  [15,  16,  17]  indicate  that  adiabatic  shear  banding  is  the  dominant  mode  of 
plugging,  erosion  and  fragmentation  in  impacts  of  4340  steel  or  depleted  uranium 
(DU)  rods  on  idled  homogeneous  armor  (RHA)  plates  over  a  range  of  velocities 
from  below  to  well  above  the  ballistic  limit.  Thus,  in  order  to  estimate  the  dis¬ 
tribution  of  fragment  sizes  and  velocities  (an  important  capability  for  terminal 
effects  code.")  it  is  essential  that  models  be  based  on  the  actual  damage  processes 
•leading  to  fragmentation. 

1.3.1  Observations  of  Shaped- Charge  Jet  Penetration  into  Steel  Ar¬ 
mor 

Shear  banding  is  always  observed  in  rolled  homogeneous  armor  (RHA)  target  ma- 
teriaJ  following  perforation  by  a  shaped-charge  jet.  Figure  8  is  an  SEM  micrograph 
showing  a  shear  band  in  a  12.7  mm-thick  RHA  plate  that  has  been  perforated  by 
a  jet  from  a  shaped  charge  warhead.  The  warhead  was  fired  at  a  standoff  of  12 
charge  diameters‘*(CD)  and  the  jet  tip  speed  was  7.7  km/s  prior  to  impact.  This 
band  has  been  exposed  by  a  diametral  section  through  the  perforation  hole.  It 
emanates  from  that  hole’s  boundary.  Three  issues  suggested  by  this  and  other 

''Throughout  this  subsection  one  charge  diameter  is  81.3  mm. 
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experiments  conducted  at  6RL  will  l>e  discussed:  structure  within  the  band;  the 
occurrence  of  a  reproducible,  characteristic  shear  band  geometry;  and  the  occur¬ 
rence  of  cracks  and/or  voids  within  the  band. 

Structure  Within  a  Shear  Band  in  Steel' Armor  At  the  magnification 
of  Figure  8,  no  structure  is  apparent  inside  the  shear  band.  This  is  in  contrast  to 
the  surrounding  material  which  still  exhibits  the  lath  structure  characteristic  of 
the  original  martensite.  (Note  however,  that  these  laths  clearly  have  been  rotated 
in  the  surrounding  material,  and  have  become  roughly  aligned  with  the  shear 
bands.)  Rogers  and  Shastry  called  white-etching  bands,  such  as  that  of  Figure  8, 
“transformed  shear  bands” ,  implying  a  phase  transformation  that  they  postulated 
to  have  occurred,  involving  austenitization  and  subsequent  quenching  [18].  They 
supported  this  hypothesis  with  measurements  of  increased  microhardness  within 
the  baind.  In  a  recent  work,  Beatty  ct  al.  [19],  presented  Selected  Area  Diffraction 
Patterns  (SADP)' within  bands  and  in  the  surrounding  martensite  in  4340  steel. 
The  diffraction  pattern  was  found  to  vary  gradually  with  distance  from  the  center 
of  the  shear  band,  with  no  clear  discontinuity  occurring  at  the  band’s  periphery. 
This  and  their  inability  to  find  remaining  austenite  within  the  bands  casts  some 
doubt  on  the  assumption  of  phase  transformation.  The  identity  of  the  structure 
within  a  shear  band  and  the  role  of  phase  transformation  in  the  genesis  of  the 
band  seem  to  be  uures*.  'ed  at  present. 

Shear  Band  Geometry  in  Steel  Armor  Our  knowledge  of  the  range 
of  geometrical  complexity  that  shear  bands  can  exhibit  is  very  limited.  Three- 
dimensional  mappings  of  shear  bands  in  target  material  would  be  useful  at  this 
time.  The  single  cross-sectional  view  of  a  shear  band  in  Figure  8  has  two  obvi¬ 
ous  geometric  features;  a  fairly  uniform  6  fim  width  £Uid  the  existence  of  little 
curvature.  These  two  observations  have  some  general  validity  in  that  they  have 
^been  found  applicable  to  other  shear  bands  emanating  from  the  target  hole  and 
exposed  by  other  radial  cuts  in  this  same  experiment.  However,  geometries  of 
shear  bands  that  formed  radially  closer  to  the  shot  line  are  difficult  to  determine 
in  such  a  ballistic  test.  The  loading  history  in  the  steel  during  its  penetration 
varied  with  position,  and  it  is  quite  possible  that  any  characteristic  features  of 
shear  band  geometry  could  have  varied  correspondingly. 

One  approach  to  the  question  of  the  sensitivity  of  shear  band  geometry  to  the 
ambient  stress  environment  during  ballistic  penetration  is  to  recover  fragments 
from  ballistic  tests.  Figures  9  and  10  2u:e  photomicrographs  which  display  a  single 
RHA  fragment  recovered  from  a  second  test  involving  perforation  of  a  12.7-mm- 
thick  plate  by  a  copper  jet.  Here  the  shaped  charge  warhead  was  fired  at  a  shorter 
standoff  of  3  CD.  The  target  plate  hole’s  fin^ll  diameter  is  14  mm,  the  same  as 
in  the  first  test  at  the  longer  standoff.  Thus  the  fragment  in  Figures  9  and  10 
indicates  the  condition  of  steel  radially  closer  to  the  shot  line  than  that  of  Figure 
8.  Figure  9  shows  the  shear  band  still  to  be  straight,  but  Figure  10  shows  its  width 
to  be  20  to  30  This  width  is  larger  than  that  of  the  sheeir  band  in  Figure 


8  by  a  factor  of  3  to  5.  Figure  11  shows  a  photomicro^aph  of  another  RHA 
fragment  from  this  test  involving  perforation  by  a  shaped-charge  fired  at  12  CD 
st^doff.  This  fragment’s  boundary  is  seen  to  be  lined  by  a  white-etdiipg  region 
of  perhaps  25  to  100  pm  width,  that  bears  little  geometrical  resemblance  to  the 
'  shear  band  of  Figure  8.  However,  the  comparison  is  marred  by  the  fact  that  the 
fragment  of  Figures  9  and  10  perforated  one  or  more  sheets  of  mild  steel  before 
coming  to  r^t  in  the  witness  pack  used  for  its  recovery.  This  points  up  the  need 
to  develop  a  non-destructive  fragment  recovery  technique  for  use  in  these  ballistic 
■  tests.  A  desirable  but  difficult- to- achie\’e  feature  of  such  a  technique  would  be 
the  preservation  of  information  regarding  each  fragment’s  original  location  on  the 
target  plate. 

A  second  approach  to  examining  the  question  of  the  dependence  of  sheau'  band 
geometry  on  original  location  in  the  target  plate  is  through  laboratorj'  experi¬ 
ments.  As  a  promising  example,  reference  [19]  documents  the  application  of  a 
split  Hopkinson  bar  technique  to  a  “hat-shaped”  specimen  of  RHA,  Labiuatory 
techniques  such  as  this  offer  a  distinct  advantage  over  baJlistic  penetration  tests  in 
that  in  the  former,  loading  conditions  within  the  specimen  can  be  relatively  well 
characterized.  However,  if  the  ultimate  goal  is  to  illuminate  phenomena  pertain¬ 
ing  to  shear  bands  that  occur  in  ballistic  penetration,  then  the  loading  conditions 
produced  in  a  laboratory  experiment  should  be  similar  to  those  occurring  at  some 
location  of  interest  in  a  target.  As  a  check,  shear  bauds  produced  in  the  laboratory 
should  be  shown  to  resemble  those  displayed  in  micrographs  from  ballistic  tests. 

Shear  Band  Association  with  Cracks  and  Voids  in  Steel  Armor  The 
ability  to  predict  information  on  fragments  produced  by  a  given  penetrator-target 
interaction  is  a  current  challenge  of  great  practical  importance.  In  this  regard, 
more  knowledge  of  the  micromechanical  processes  leading  to  fragmentation  is 
needed.  Micrographs  such  as  that  of  Figure  12  seem  to  implicate  shear  bands. 
This  shows  a  ^hear  band  in  the  same  target  plate  considered  in  Figure  8.  A 
crack  clearly  extends  along  part  of  the  shecu  band’s  length.  Figure  13  shows 
the  same  shear  band  and  crack  at  smaller  magnification.  The  crack  is  seen  to 
bepn  at  the  hole  boundary;  the  process  of  fragment  formation  has  apparently 
been  caught  in  a  state  of  partial  completion.  The  sequence  of  events  linking  shear 
band  localization  to  macroscopic  cracking  is  still  unclear.  Perhaps  the  increased 
microhardness  within  the  band  provides  a  clu 

Spheroidal  voids  are  shown  in  the  same  target  plate  in  Figure  14.  These 
have  coalesced  to  form  a  crack,  along  which  the  martensitic  laths  have  become 
aligned.  Such  voids  may  constitute  a  second  micromechanical  mechanism  for 
fragmentation.  However,  it  is  not  yet  clear  to  .what  extent  they  axe  indeed  a 
separate  mechanism,  independent  of  shear  banding 


1.3.2  Comparison  of  Uranium  Alloys  and  Tungsten  Composities  as 
Long  Rod  Penetrator  Materials 

Tlie  importance  of  both  plastic  inhomogeneities  and  material  anisotropies  to  the 
process  of  penetration  is  also  illustrated  by  a  comparison  of  the  observed  phe¬ 
nomenology  amd  resulting  ballistic  performance  of  uranium  alloy,  tungsten  heavy 
alloy  composite,  and  tungstten  single  crystal  penetrators. 

Depleted  uranium  (DU)  alloys  and  tungsten  heavy  alloys  (WHA)  offer  both 
the  high  materiaJ  densities  (for  penetration  performance)  and  useful  mechanical 
properties  (needed  for  the  integration  of  the  penetrator  core  in  engineered  ki- 
netic  energj'  (KE)  long  rod  projectiles  and  cartridge  assemblies).  MetaJlurgically, 
however,  DU  amd  WHA  are  very  different  kinds  of  materiads. 

The  DU  materials  aire  true  aJloys.  The  standaird  DU-3/4Ti  alloy  cam  be  caist 
and  heat-treated  in  much  the  same  way  as  iron-carbon  (steel)  aJloys,  Eckelmeyer 
(20).  In  contrast.  Tungsten  is  a  refractory  metal  and  pure  tungsten  products 
generally  have  poor  ductilities  amd  toughnesses,  and  require  very  high  processing 
temperatures.  Therefore  the  WHAs  used  for  KE  penetrator  applications  are  actu- 
aJly  composites  of  tungsten  particles^  in  a  ductile  metaJ  matrix  (usually  consisting 
primarily  of  nickel' and  iron),  and  usually  produced  hy  liquid-phaise  sintering  tech¬ 
niques.  The  density  of  the  WHA  composite  depends  on  its  tungsten  content,  aind 
the  mechanical  properties  depend  on  both  its  tungsten  content  (percentage  of 
ductile  metal  matrix)  and  on  its  processing  history,  Gurwell  et  al.  [21].  Because 
DU  amd  WHA  aire  so  very  different  metallurgically,  they  also  differ  significantly 
in  many  fundamental  mechanicail,  thermal,  etc.  properties. 

The  ballistic  performamces  of  DU  amd  WHA  projectiles  aJso  differ  by  a  small, 
but  significamt,  aunount.  Penetrators  of  the  standard  DU-3/4Ti  alloy  consistently 
require  lower  impact  velocities  (by  approximately  100  m/s)®  to  defeat  thick  steel 
amd  steel/ceramic  launinate  targets  than  WHA  projectiles  of  equal  mass  and  iden- 
ticad  geometry  (i.e.  equal  density)  Magness  and  Fairramd  [22]. 

Exauninations  of  residual  DU  and  WHA  penetrators  embedded  in  steel  armor 
(Figures  15  amd  16,  respectively)  reveail  that  the  tv/o  materiads  display  very  dif¬ 
ferent  deformation  behaviors  during  the  penetration  process.  In  the  caise  of  the 
WHA  penetrator  (Figure  l6),  the  relatively  smooth  erosion  amd  back-fextrusion 
of  the  material  at  the  penetrator’s  head  is  quaditatively  consistent  with  the  fluid- 
dynamic  idealization  of  penetration  depicted  in  Figure  1.  The  WHA  material 
undergoes  a  great  deal  of  plastic  deformation  as  it  approaches  and  is  inverted  at 
the  penetrator-target  interface.  A  laurge  mushroomed  head  is  found  on  the  residuad 
penetrator  embedded  in  the  tairget  (Figure  16)  amd  is  routinely  observed  on  resid¬ 
ual  WHA  penetrators  aifter  exiting  finite  steel  targets  (Figure  17).  The  evolution 
of  the  WHA  microstructure  in  Figure  16,  from  the  originadly  sphericad  tungsten 

®This  relatively  small  velocity  difference  translates  into  a  significant  difference  («  2  km)  in 
^the  effective  range  of  such  projectiles  by  virtue  of  air  drag  effects. 


pa? tides  to  highly  elongated  stringers,  provides  further  evidence  of  the  large  plas¬ 
tic  strains  that  occur  as  the  material  is  inverted.  On  closer  examination,  of  course, 
the  flow  of  the  penetrator  materid  does  differ  from  fluid  behavior.  Because  the 
WHA  is  a  solid,  small  elastic  strains  and  some  small  plastic  strains  manage  to 
propagate  away  from  the  interface.  The  “flaring”  of  the  residual  rod  as  it  feeds 
into  the  mushroomed  head  is  evidence  of  this.  Another  significant  departure  from 
the  fluid-dynamic  ideal  is  the  development  of  plastic  localizations,  or  small  re¬ 
pons  of  larger  than  average  plastic  strains,  near  the  periphery  of  the  mushroom. 
The  eventual  discard  of  material  from  the  periphery  of  the  mushroomed  head 
ultimately  occurs  along  some  of  these  localizations. 

The  deformation  of  the  DU-3/4Ti  penetrator  has  a  very  different  appearance 
(Figure  15).  Localizations  of  the  plastic  deformation  develop  very  quickly  along 
planes  of  highest  local  shear  stress  and  at  relatively  small  values  of  strain.  The 
subsequent  discard  of  the  back-extruding  material  along  these  shear  planes  occurs 
before  a  large  mushroomed  head  can  develop  on  the  penetrator.  -^s  a  result, 
the  residual  DU  penetrators  observed  exiting  the  rear  face  of  finite  st^k-t^gets 
lack  the  large  mushroomed  head  of  the  WHA  penetrators.  Instead,  they  display 
a  pointed  nose  (Figure  18),  shaped  by  the  chiseling  action  of  the  local  shear 
f^Hures  occurring  during  the  penetration  process.  It  is  the  ability  of  the  DU  alloy 
penetrators  prevent  the  build-up  of  a  large  mushroomed  head  which  is  ultimately 
responsible  for  their  superior  ballistic  performance. 

A  comparison  of  the  penetration  cavities  formed  by  DU  and  WHA  projectiles 
penetrating  into  thick  steel  blocks  is  shown  in  Figure  19.  The  hole  profiles  in 
Figure  19  reveal  that,  during  much  of  the  penetration  process,  the  DU  penetrator 
creates  a  smaller  diauneter  penetration  tunnel.  Therefore,  less  of  the  kinetic  energy 
of  the  DU  projectile  is  expended  laterally,  to  displace  a  l^ger  diameter  hole  into 
the  target,  and  more  of  its  kinetic  energy  is  available  to  penetrate  to  a  greater 
depth  in  the  thick  target.  Against  lesoer  thicknesses  of  armor  plate,  which  both 
the  DU  and  WHA  projectiles  are  able  to  perforate,  the  DtJ  penetrator  does  so 
more  efficiently  and  therefore  has  a  lower  ballistic  limit.® 

The  distinct  deformation  behaviors  exhibited  by  the  DU  and  WHA  penetrators 
are  a  result  of  differences  in  their  ffiermomechanical  properties.  Under  the  very 
high  strain  rate  loading  conditions  of  the  penetration  process,  there  is  essentially 
no  time  for  the  macroscopic  transport  of  the  heat  being  generated  by  the  defor¬ 
mation  of  the  penetrator.  The  thermal-softening  induced  by  this  adiabatic  heat¬ 
ing  overcomes  the  usual  strain-hardening  and  strain-rate  hardening  much  more 
quickly  in  the  DU  than  in  the  WHA.  The  result  is  a  net  strain-softening  behavior 
in  the  DU,  under  which  small  perturbations  in  the  plastic  strain  or  temperature 
fields  can  quickly  grow  into  the  plastic  localization.s  (adiabatic  shear  bands)  ob- 

*For  a  given  penetrator  and  target,  the  ballistic  limit  is  defined  to  be  the  lowest  penetrator 
velocity  which  will  produce  perforation. 


served  in  the  DU  residual  penetrators.  As  evidenced  by  the  smoothly  deformed 
microstructures  and  the  large  mushroomed  heads  which  develop  on  the  WHA 
penetrators,  the  plastic  deformation  is  much  more  stable  in  this  material. 

A  related,  but  distinctly  difiFerent,  illustration  of  the  role  of  the  penetrator 
material’s  behavior  on  penetration  performsmce  is  the  comparison  of  the  behavior 
and  performcince  of  different  orientations  of  tungsten  monocryst2j  penetrators  and 
polycrystalline  pure  tungsten  penetrators.  Due  to  the  thermomechanical  proper¬ 
ties  of  pure  tungsten  (high  strain  rate  sensitivity,  relatively  low  rate  of  thermal 
softening,  etc.),  the  plastic  deformation  of  the  tungsten  polycrystal  penetrator 
,is  extremely  stable  (Figure  20).  In  fact,  the  “mushrooming”  of  the  tungsten 
polycrystal  is  more  extensive  than  that  of  the  WHA  resulting  in  poorer  ballistic 
performance  than  a  WHA  penetrator  of  equivaJent  mass  and  geometry,  Mag- 
ness  [22].  However,  ballistic  tests  show  that  the  performance  of  monocrystalline 
tungsten  penetrators  (Bruchey,  Horwath  and  Kingman  [23])  vary  with  their  crys¬ 
tallographic  orientations.  In  one  case,  with  the  axis  of  the  penetrator  parallel  to 
the  <110>  direction  in  the  lattice,  the  mushrooming  can  be  even  more  exagger¬ 
ated  than  in  the  pure  tungsten  polycrystal  (Figure  21).  As  a  result,  the  diameter 
of  the  penetration  cavity  is  considerably  larger  and  the  projectile  achieves  very' 
low  depths  of  penietr''tion.  In  contrast,  if  the  axis  of  the  rod  is  parallel  to  the 
<100>  axis  of  the  crystal,  cleavage  failures  aippear  in  the  head  of  the  penetrator 
after  only  small  plastic  strains  (Figure  22)  and  the  deformed  material  is  more 
quickly  discarded.  The  <100>  oriented  crystal  penetrator  creates  a  much  smaller 
diameter  penetration  tunnel  and  delivers  penetration  performance  approximating 
that  of  DU.  In  this  case,  anisotropies  in  the  strengths  and  flow  behaviors,  rather 
than  plastic  flow  localization,  lead  to  the  very  different  modes  of  deformation  and 
ballistic  performances  exhibited  by  the  monocrystaJline  tungsten  penetrators. 

The  differing  behaviors  and  bcdlistic  perforniances  of  DU,  WHA,  and  single 
crystal  tungken  projectiles  serve  as  excellent  tests  of  the  ability  of  ballistic  impact 
simulation  codes  to  Jiccurately  represent  the  response  of  penetrator  and  armor  ma¬ 
terials  during  the  penetration  event.  Conversely,  the  incorporation  into  simulation 
codes  of  state-of-the-art  materieJ  models  and  algorithms  should  become. a  central 
goal  for  the  code  development  community. 

1.3.3  Conclusions 

The  fundamental  reason  for  non-hydrodynamic  behavior  in  the  penetration  of 
metallic  targets  is  the  complex  internal  microstructure  of  melals  2uid  its  mani¬ 
festations  as  macroscopic  strength  and  localized  deformation  modes  which  lead 
to  material  failure.  The  occurrence  of  ductile  faiilure  modes  depends  strongly  on 
several  aspects  of  the  ballistic  impact. 

.  „  .  .  .  *  . , 

1.  The  macroscopic  theimo-visco-plastic  (TVP)  response  of  the  material. 


2.  Features  of  the  material  microstnicture  which  can  serve  as  nucleation  sites 
for  shear  brad  or  voids. 

3.  Features  of  the  macroscopic  loading  (e.g.  hydrodynamic  pressure/tension, 
stress  concentrations,  etc.)  which  can  encourage  or  supress  localized  defor¬ 
mation. 

The  macroscopic  TV?  response  of  high  strength  steels  is  quite  favorable  for  shear 
band  formation  and  in  this  context  some  aspects  of  band  morphology  and  met¬ 
allurgy  were  discussed  in  section  1.3.1.  Most  notable  is  the  severity  of  the  lo¬ 
calization  in  a  shear  band  and  the  fact  that  bands  may  serve  as  nucleators  for 
other  ftulure  mechanisms  e.g  cracks.  Lastly,  section  1 .3.2  emphasized  the  degree 
to  which  a  material’s  shear  band  susceptibility  is  affected  by  its  TVP  response. 
Moreover,  the  relative  shear  band  susceptibility  of  different  materials  produced 
markedly  different  penetration  performance. 

2  Terminal  Effects  Codes  (K.  Kimsey,  G.  Randers- 
■  Pehrson) 

Computer  codes  used  to  study  the  armor  penetration  process  are  usually  known  as 
“hydrocodes”,  or  more  accurately,  “wave  propagation  codes’".  Wave  propagation 
codes  fall  into  two  classes:  Etilerian  .and  Lagrabgian  (a  third  class,  Arbitrary 
Lagrangi^m  Eulerian,  is  currently  emerging).  Eulerian  codes  used  at  BRL  include 
HULL  [24],  MESA  [25],  and  CTH  [26].  Lagrangian  codes  used  at  BRL  include 
DYNA2D/3D  [27]  and  EPIC2D/3D  [28]  and  some  of  their  derivatives.  CTH  and 
MESA  are  now  being  evaluated  for  modeling  arir  or /anti-armor  problems. 

In  either  approach,  the  problem  is  divided  into  a  large  number  of  small  compu¬ 
tational  .cells  or  elements.  The  solution  to  the  prablem  proceeds  over  small  time 
steps,  typically  a  tenth  or  hundredth  of  a  microsecond.  At  each  time  step  the 
stress  state  within  each  computationcJ  cell  is  determined  from  the  deformation 
of  the  cell  together  with  the  equation  of  state  (pr  assure- volume  relationship)  and 
constitutive  model  (describing  the  thermo- visco-j  lastic  response)  of  the  material 
that  it  contains.  From  the  stress  state,  the  equa  ion  of  motion  is  used  to  deter¬ 
mine  the  defonnation  to  be  used  as  input  to  the  next  time  step.  The  Lagrangian 
and  Eulerian  approaches  differ  in  the  treatment  of  material  motion.  In  the  La¬ 
grangian  codes,  the  cell’s  mass  is  assigned  to  the  nodes  at  the  cell  comers,  aind  the 
nodes  move  through  space  according  to  the  equation  of  motion.  The  Lagrangiaua 
cell  becomes  deformed  and  moves  along  with  its  comer  nodes.  In  the  Eulerian 
approach,  the  cell  does  noLmove,  but  material  from  reighboring  cells  flows  in  or 
out  according  to  the  equation  of  motion.  All  of  the  Eulerian  codes  use  a  two-step 
solution  scheme.  The  first  step  is  a  Lagramgian  si  ep  in  which  the  cells  distort  to 
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follow  material  motion.  The  second  step  is  a  remesh  step  where  distorted  cells  are 
mapped  back,  i.e.  advected,  to  the  Eulerian  mesh. 

Because  of  the  severe  deformations  which  occur  in  the  penetration  process, 
Lagrangian  codes  generally  cannot  treat  such  problems'  unless  special  measures 
are  taken.  In  EPIC  and  DYNA,  “eroding  sliding  interfaces”  are  used,  wherein 
severely  deformed  cells  are  removed  from  the  computation  ax:cording  to  criteria 
such  as  accumulated  plastic  strain.  Such  special  treatments  extend  the  utility  of 
Lagrimgian  codes  for  modeling  deep  penetration  problems.  Eulerian  codes  can 
model  severe  deformation  more  readily,  without  resorting  to  special  treatments, 
but  they  too  have  disadvantages,  for  example  diflBculties  in  advecting  material 
stress  histories.  Eulerian  codes  generally  require  much  larger  amounts  of  computer 
memory  and  processor  time. 

Most  armor/anti-£irmor  problems  are  inherently  three-dimensiond.  Some  two- 
dimensional  problems  such  as  normd  impact  and  penetration  of  semi-infinite  ar¬ 
mor  are  useful,  however,  mainly  for  validation  of  odes  Following  are  brief 
descriptions  of  several  three-dimensional  wav.  propagation  codes  in  current  use 
for  penetration  mechanics  studies  at  iiRL. 

HULL  is  a  family  of  computer  programs  for  solving  continuum  dynamics  prob¬ 
lems.  The  code  solves  the  finite  clifferciice  analoiis  of  the  governing  partid  differen- 
tid  equations  for  two-  and  three-dimensiond  Eulerian  and  Lagrangian  formulated 
problems  as  well  as  coupled  Eulerian/Lagrangian  solutions.  Materid  advection  is 
based  on  a  fixst-order  donor-cell  dgorithm  with  a  heuristic  multimaterid  diffusion 
limiter  to  preserve  material  interfaces.  Most  solids  in  HULL  axe  modeled  with  the 
Mie- Gruneisen  equation  of  state.  Concrete,  geologicd  materids  and  explosive 
products  equations  of  state  are  dso  avdlable.  Materid  fdlure  models  include 
rvaximum  principd  stress,  maximum  principd  strdn,and  the  Hancock-Mackenzie 
t  'iaxid  fdlure  model.  Plasticity  models  include  elastic-perfectly  plastic  with  pro- 
v’siotiSi  for  work  hardening  and  thermal  softening. 

MESA  is  a  new  3D  Eulerian  code  which  is  being  developed  for  armor/anti- 
armor  work  by  a  team  at  Los  Alarhos  Nationd  Laboratory.  This  code  was  written 
to  take  advantage  of  new  high-order  accurate  advection  dgorithms  and  interface 
reconstruction  methods.  Plasticity  models  in  MESA  include  the  Johnson-Cook 
and  Steinberg- Cochran- Guinan- Lund  models,  and  various  simpler  ones.  The  duc¬ 
tile  fracture  model  TEPLA-F  [29]  has  dso  been  incorporated  in  MESA. 

CTH  is  Sandia’s  new  Eulerian  code.  Like  MESA,  it  incorporates  new  high- 
order  accurate  dgorithms.  It  h^  its  roots  in  the  previous  Sandia  code,  CSQ. 
CTH  solves  the  finite  difference  analogs  of  the  governing  conservation  equations 
using  an  Eulerian  solution  scheme.  CTH  models  solid  dynamics  problems  involv¬ 
ing  shock  wave  propagation,  multiple  materids  and  large  deformations  in  one, 
two  and  three  dimensions.  CTH  models  elastic-plastic  behavior,  high  explosive 
detonation,  fracture  and  motion  of  fragments  smdler  than  a  computatioad  cell. 
Plasticity  models  include  elastic-perfectly  plastic  with  thermd  softening,  Johnson- 
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OK>k,  Zerilli-Annstrong  and  Steinberg-Guinan-Lund.  The  Johnson-Holmquist 
[30]  brittle  material  model  has  also  been  incorporated  in  CTH.  Both  analytic^ 
(Mie-Gruneisen)  and  tabular  equations'  of  state  are  available  for  modeling  the 
hydrodynamic  behavior  of  materials.  Three  types  of  high  explosi^  bum  model 
are  avmlable;  programmed  bum,  CJ  volume  bum  and  history  variable  reactive 
bum.  The  Jones- Wilkins-Lee  equation  of  state  is  available  for  modeling  detona¬ 
tion  products.  Fracture  initiation  is  based  on  maximum  pressure  and  principal 
stress. 

DYNA3D  is  a  Lagrangian  explicit  finite  element  code  for  analyzing  the  tran¬ 
sient  dynamic  response  of  three-dimensional  solids  and  structures.  A  product  of 
the  Lawrence  Livermore  National  Laboratory,  it  has  a  wide  variety  of  element 
types  and  material  models  and  sophisticated  treatment  of  contact/impact  inter¬ 
faces.  A  version  of  the  code  supported  by  Dr.  John  Hall quist,  the  original  author  of 
the  code  and  now  an  independent  consultant,  incorporates  an  automatic  eroding 
interface  algorithm.  Recently,  researchers  at  LLNL  have  implemented  a  slideline 
adaptive  node  definition  (SAND),  i.e..  eroding  slidelines,  algorithm  into  LLNL’s 
DYNA3D  code. 

EPIC  is  a  Lagrangian  explicit  finite  element  code  written  by  Dr.  Gordon 
Johnson  of  Alliant  Tech  Systems  (formerly  a  division  of  Honeywell).  It  incorpo¬ 
rates  an  eroding  surface  capability,  and  uses  primarily  Johnson’s  material  models 
(Johnson-Cook  for  metals,  Johnson-Holmquist  for  ceramics,  and  the  JWL  equa¬ 
tion  of  state  for  detonation  products).  The  code  uses  triangular  or  tetrahedral 
finite  elements  rather  than  quadrilateral  or  brick-shaped  ones,  and  consequently 
can  frequently  run  problems  with  large  deformations  more  readily  than  other  La¬ 
grangian  codes. 

ALE  (Arbitrary  Lagrangian  Eulcrian)  codes  under  development  include  the  2D 
CALE  code,  written  in  the  C  programing  language  at  Livermore,  and  a  3D  mixed- 
zone  code  being  developed  by  Livermore.  The  latter  code  is  written  in  Fortran 
and  has  its  roots  in  the  DYNA3D  code.  Like  DYNA3D,  it  has  a  wide  variety 
of  plasticity  models,  including  Johnson-Cook,  Steinberg-  Cochran- Guinan-Lund, 
power  law  hardening,  crushable  foam,  geologic  models,  and  various  orthotropic 
models.  Sandia  is  also  working  on  a  3D  ALE  code  called  RHALE. 

In  this  article,  we  focus  on  the  material  modeling  and  equation  of  state  mod¬ 
eling  in  present  day  ternndnal  ballistics  codes.  The  predictive  capability  of  current 
wave  propagation  codes  for  impact  studies  is  dependent  on  the  materid  model 
and  properties  that  are  used  in  the  Simulations.  Metals,  until  recently,  were  the 
materials  of  primary  interest.  It  is  common  in  existing  wave  propagation  codes 
for  impact  studies  to  divide  the  deformation  behavior  of  metals  into  volumetric 
and  deviatoric  parts  as 

Cij  =  -p6,j  -I-  ,  (2) 

in  which  <r  is  the  stress  tensor,  s  is  its  deviatoric  part  and  p  is  the  pressure.  The 
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volumetric  behavior  is  obtained  from  an  equation  of  state  which  describes  the  in¬ 
terdependence  of  thermodynamic  properties,  such  as  pressure,  energy,  density  and 
temperature.  For  impacts  where  the  striking  velocity  is  below  the  sonic  velocity, 
some  variant  of  the  Mie-Gruneisen  equation  of  state  is  commonly  used.  Beyond 
this  range  (hypervelocity  impact)  the  Tillotson  equation  of  state  is  favored.  Some 
codes  also  use  tabular  equations  of  state  to  represent  equation  of  state  surfaces 
of  vbitrary  complexity  which  are  valid  for  a  wide  range  of  densities  and  temper¬ 
atures.  Examples  of  tabular  equations  of  state  include  SESAME  (Bennett  et  al. 
[31],  Kerley  [32])  and  ANEQS  (Thompson,  [33]).  Explosive  detonation  products 
are  most  commonly  represented  by  the  well-known  JWL  amalytical  equation  of 
state.  A  desirable  feature  of  the  analytical  equation-of-state  models  is  that  they 
are  easy  to  use  and,  for  most  metals,  model  parameters  can  be  obtained  from 
the  literature  (van  Thiel  [34];  Marsh  [35];  Kohn  [36];  Steinberg  [37];  Dobratz  and 
Crawford  [38]). 

An  incremental  elastic-pltistic  formulation  is  used  to  describe  the  deviatoric 
response  of  metals  in^^esent  finite  difference  and  finite  element  codes.  Most 
follow  the  description  first  given  by  Wilkins  [39,  40].  The  plasticity  descriptions 
are  usually  based  on  the  assumed  decomposition  of  the  strain  rate  tensor  into 
elastic  2md  plastic  parts 

=  +  (3) 

together  with  incompressibility  of  the  plastic  part 

£<5  =  0.  (4) 

The  von  Mises  yield  criterion  is  commonly  used  to  determine  if  a  material  element 
is  imderjoing  plastic  flow. 

Traditionally,  the  plasticity  model  -  in  wave  propagation  codes  for  impact  stud¬ 
ies  have  been  relatively  simple,  of  the  elastic,  perfectly-plastic  type,  and  have  been 
implemented  using  the  radied  return  algorithm  formulated  by  Wilkins  [39,  40]. 
More  recently  viscoplastic  models  have  been  implemented  in  wave  propagation 
.'odes  to  represent  the  dependence  of  the  flow  stress  on  strain,  strain  rate,  pres¬ 
sure,  temperature  or  some  combination  thereof.  Viscoplastic  models  for  metals 
include  those  of  Johnson  and  Cook  [41]  and  Zerilli  and  Armstrong  [42]. 

The  Johnson- Cook  model  for  the  von  Mises  equivalent  flow  stress  is  expressed 

r  . 

(T,  =  [A.bB(e»’)"]  l.fCln^^^^j[l-(r-)"'],  (5) 

where  ctj  is  the  flow  stress,  (P  is  the  equivalent  plastic  strain,  fP  is  the  rate  of 
change  of  the  latter,  and  T*  is  the  homologous  temperature;  A,  B,  n,  C  and  m 
are  constants.  The  expression  in  the  first  pair  of  brackets  represents  the  depen¬ 
dence  of  the  flow  stress  on  (von  Mises)  equivalent  strain.  The  expressions  in  the 
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a^^iid  and  tliird  brackets  reprraent  the  dependence^t  strain  rate  and  tempera- 
tu^,  ^pectivdy.  The  model  is  based  on  a  numerical  fit  to  torsion  and  tension 
tot  data  at  ambient  and  elevated  temperatures.  Use  of  the  model  outside  the 
ran^  of  the  test  data  may  not  be  warranted.  Model  parameters  for  a  number 
oi  metals  am  given  by  Johnson  and  Cook  [41]  and  Johnson  and  Holmquist  [43]. 
The  JohusomCook  mudel  has  been  implemented  in  the  EPIC,  DYNA,  CTH  and 
MESA  a>des. 

Figure  23  compares  computed  semi-infinite  penetration  depths  with  quarter- 
scale  test  data,  Magness  [44],  for  65  gram,  L/D  10,  depleted-urcinium  alloy  (U- 
3/4Ti)  and  tungsten  heavy  alloy  (WHA)  rods  impacting  semi-infinite  RHA  targets 
at  striking  velocities  between  800  and  1800  m/s.  The  solid  s)nmbols  represent 
CTH  results  using  the  Johnsdn-Cook  model  to  describe  the  flow  stress  for  both 
rod  Md  target  materials.  The  CTH  results  compare  well  with  the  test  data 
and  iJso  reflect  the  observed  gap  in  performance  between  these  two  materials. 
Magness  ?uid  Farrand  [22]  conducted  ballistic  tests  and  microstructural  analyses 
of  the  performance  of  uranium  and  tungsten  heavy  alloys  and  identified  the  strain¬ 
softening  behavior  developed  by  the  uranium  alloys  under  high  rate  loading  and 
the  resulting  thermo-mechanical  (adiabatic  shear)  instability,  as  the  mechanism 
responsible  for  uranium’s  superior  performance.  It  is  interesting  to  note  that 
'while  the  Johnson-Cook  model  does  not  explicitly  model  adiabatic  she<ir  bands,  it 
appears  to  capture  some  aspects  of  the  “global”  behavior  of  the  two  alloys  under 
impact  loading  (for  the  striking  velocities  studied). 

.  Zerilli  and  Armstrong  [42]  have  shown  an  improved  description  of  copper  and 
iron  cylinder  impact  (Taylor)  test  results  using  their  dislocation-mechanics-based 
constitutive  relations  for  body  centered  cubic  (bcc)  and  face  centered  cubic  (fee) 
materials.  The  Zerilli-Armstrong  model  is  based  on  the  observation  that  each 
material  structure  type  (fee,  bcc,  hep)  will  have  its  own  constitutive  behavior, 
dependent  on  the  dislocation  characteristics  for  a  particular  structure.  Thus, 
separate  relations  are  required  to  describe  flow  stress  for  bcc  and  fee  metals.  For 
OFHC  copper  (fee),  the  flow  stress  is  given  by 

<r  =  AtTa'  +  C2e*/*exp[(-‘C3  -I-  In  e)r]  -f  ,  (6) 

whereas  for  Armco  iron  (bcc)  the  flow  stress  is  given  by 

=  Ciexp[(-<^-|-c<ine)T]-f  C5e"-t-Arf-’/’;  (7) 

the  constitutive  variables  have  the  same  interpretation  as  in  (5)  except  that  t  is 
the  grain  size.  Values  for  the  eight  material  constants:  Cj, . . .  ,05,  Affg',  k,  n,  for 
OFHC  copper  and  Armco  iron  are  reported  by  Zerilli  and  Armstrong  [42];  this 
model  is  currently  available  in  EPIC,  CTH  and  MESA. 

The  Mechanical  Threshold  Stress,  Follansbee  [45],  and  the  Steinberg- Cochran- 
Guin^-Lund,  Steinberg  et  al.  [46,  47],  material  models  relate  the  macroscopic 
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stress-strain  behavior  to  Tiicroscopic  properties.  Material  parameters  have  been 
published  for  the  latter  model  for  seven  metals:  beryllium,  molybdenum,  niobium, 
niobium  alloy  fansteel  85,  tantalum,  tungsten,  vanadium,  and  for  some  other 
materials  (Steinberg  [37]).  Both  models  have  been  implemented  in  MESA. 

Computational  failure  models  for  impact  loading  situations  are  discussed  in 
review  articles  by  Seaman  [48]  and  Zukas  [49]  and  the  NMAB  Committee  re¬ 
port  on  Materials  Response  to  Ultra-High  loading  Rates  [1].  Briefly,  empirical 
failure  models  of  varying  degrees  of  complexity  exist.  Some  have  been  applied  sue-, 
cessfully  in  high  velocity  impact  simulations.  For  the  most  part,  though,  failure 
criteria  and  models  are  of  an  ad-hoc  nature,  lacking  a  micromechanical  basis  to 
comprehensively  treat  problems  involving  both  brittle  and  ductile  failure  modes. 
Different  criteria  apply  for  different  impact  scenarios  and  there  are  at  present  no 
guidelines  for  analysts  wishing  to  select  a  failure  criterion  appropriate  to  varying 
conditions. 

The  simplest  of  the  failure  initiation  criteria  are  based  on  instantaneous  values 
of  a  field  variable,  such  as  pressure,  stress,  strain,  plastic  work,  or  some  combina¬ 
tion  thereof.  Once  the  criterion  has  been  satisfied  at  a  given  location,  failure  is 
considered  to  occur  instantaneously.  The  post-failure  behavior  of  materials  can 
be  described  in  a  number  of  ways^  The  failed  material  may  be  removed  entirely 
from  the  calculation  or  be  described  by  a  modified  constitutive  function. 

Time-dependent  initiation  criteria  represent  the  next  level  of  sophistication 
and  have  been  successfully  applied  in  several  situations.  One  of  the  earliest  tiihe- 
depjendent  models  is  attributed  to  Tuler  and  Butcher  [50].  Failure  is  assumed 
to  occur  instantaneously  only  after  the  stress  has  remained  above  a  threshold 
stre  ss  level  for  some  length  of  time,  below  which  no  significant  damage  will  occur 
regardless  of  duration. 

Criteria  in  which  damage  accumulation  is  a  function  of  field  variables  have  been 
dev;loped  by  Davison  [51]  and  Johnson  [52]  In  these  models  the  rate  of  damage 
accumulation  is  taken  to  be  a  function  of  strain,  temperature,  pressure,  strain 
rat<  and  current  damage  level.  Damage  levels  are  used  to  degrade  a  materials’ 
cha -acteristics  and  sometimes  trigger  instantaneous  failure. 

flesecirchers  at  SRI  International  (Seaman  and  Shockey  [53],  Se^tin  et  al.  [54] 
and  Erlich  [55])  have  been  developing  “microstatistical”  failure  models.  Ductile 
failure  damage  is  initiated  when  the  average  stress  exceeds  a  tensile- pressure  cri¬ 
terion.  Brittle  fracture  is  initiated  when  the  maximum  norm2J  stress  exceeds  a 
tens  ile  threshold.  Shear  banding  begins  when  the  maximum  plastic  shear  strain 
exc<eds  a  criticcil  value.  After  initiation,  voids,  cracks,  and  shear  bands  nucleate 
and  grow  according  to  experimentally  determined  rate  equations.  These  failure 
moc  els  have  not  been  implemented  in  many  production  codes  due  the  limited 
number  of  materials  which  have  been  chara'-terized  as  well  as  the  requirement  to 
conduct  extensive  materials  tests  in  order  to  obtain  model  parameters  for  new 
materials. 


No  ddinite  failure  model  has  yet  emerged,  nor  has  a  library  of  data  to  drive 
the  various  models  ior  materials  of  practical  interest.  Therefore,  for  want  of 
an  idequate  data  1 4se  for  advanced  failure  models,  most  high  velocity  impact 
simulations  are  perfonned  with  the  simplest  fsulure  models.  In  many  cases,  simple 
f^lu]%  models  produce  results  in  agreement  with  experiments. 

Brittle  material  modeling  is  a  topic  currently  of  great  interest  in  the  compu¬ 
tational  penetration  mechcmics  community,  primarily  because  of  the  interest  in 
ceramics  as  candidate  armor  materials.  Despite  the  recent  surge  in  activity,  there 
does  not  yet  appear  to  be  common  agreement  in  the  developer’s  community  on 
the  correct  modeling  approach.  Nonetheless,  the  first  models  intended  to  address 
the  behavior  of  ceramic  materials  under  ballistic  impact  are  becoming  avrilable 
to  the  general  user  community. 

Probably  the  most  notable  of  these  new  models  is  the  Johnson- Holmquist 
model  |30],  which  has  been  implemented  in  the  1990  version  of  the  EPIC  code, 
Johnson  and  Stryk  f28],  The  salient  features  of  this  model  include  Mohr-Coulomb- 
■  like  behavior  of  the  yield  surface  (i.e.  yield  strength  proportional  to  pressure), 
and  an  equation  controlling  the  evolution  of  a  scalar  “damage”  parameter.  When 
the  damage  parameter  reaches  a  critical  value,  failure  is  modeled  by  imposing  a 
pressure  jump  intended  to  account  for  shear- induced  dilatancy  believed  to  occur 
in. the  comminuted  ceramic.  However,  neither  the  magnitude  nor  the  onset  of  the 
pressure  jump  is  a  function  of  material  shearing  rate,  as  would  be  indicated  by 
granular 'flow  theory. 

Each  of  the  brittle  material  models  which  axe  beginning  to  surface  in  the  com¬ 
putational  penetration  mechanics  community  address  issues  of  ceramic  response 
deemed  important  by  their  respective  developers.  The  features  of  ceramic  re¬ 
sponse  which  are  being  addressed  , in  various  forms  by  various  developers  include: 
volumetric  bulking  upon  failure;  shear  induced  dilatancy;  frictional  effects  in  failed 
ceramic  rubble;  modulus  degradation;  Mohr- Coulomb- like  failure  surfaces  in  ce¬ 
ramics;  microscopic  cracking  and  void  formation;  2uid  “damage”  evolution.  It 
remains  to  be  seen  which  of  these  model  features  aie  necessary  to  model  the  bal¬ 
listic  behavior  of  ceramic  materials.  Toward  this  end  Kipp  and  Grady  of  Sandia 
are  conducting  experiments  (mostly  one-dimensional  plate  impacts)  to  measure 
the  dynamic  response  of  ceramics. 

3  Finite  Strain  Plasticity  (N.  Huffington) 

3.1  State  of  Modeling  in  Current  Terminal  Ballistics  Codes 

It  should  be  understood  that  the  theory  of  plasticity  has  not  achieved  a  state  of 
uniform  acceptance  comparable  to  the  theory  of  elasticity,  even  for  relatively  small 
strains.  Even  such  basic  concepts  as  the  plastic  potential  and  use  of  a.ssociated 
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vs.  non-associated  flow  rules  are  still  unsettled.  For  the  historic2d  development  of 
plasticity  theory  one  may  refer  to  Hill  [56]  for  a  condensation  of  results  up  to  1950 
tmd  to  Malvern  [57]  for  publications  up  to  1969.  Subsequent  publications  (very 
numerous)  are  scattered  throughout  the  research  literature.  Since  large  plastic 
strains  are  possible  only  in  materials  which  possess  considerable  ductility,  the 
discussion  which  follows  is  essenticdly  concerned  with  various  aspects  of  constitu¬ 
tive  modeling  required  for  representation  of  behavior  of  ductile  metals  (rubberlike 
materials  are  usually  treated  using  specialized  nonlinear  elastic  constitutive  mod¬ 
els).  No  attempt  will  be  made  to  itemize  which  of  these  desirable  features  are 
currently  available  in  specific  codes;  this  information  is  usually  available  in  the 
documentation  for  each  code. 

Work  Hardening  Provision  for  work  hardening  may  vary  from  none  (per¬ 
fectly  plastic  material)  to  linecir  hardening  (one  hardening  modulus)  to  nonlinear 
hardening:  either  a  particular  functional  dependence  of  stress  on  strain  is  specified 
or  tabular  data  on  stress-strain  pairs  is  used  in  conjunction  with  an  interpolation 
routine.  Linear  hardening  is  usually  satisfactory  only  over  a  limited  strain  range 
where  the  maximum  strain  is  known  in  advance.  For  materials  such  as  copper  a 
nonlinear  representation  is  essential  since  no  portion  of  the  stress-strain  curve  is 
linear.  For  other  than  uniaxial  loading  it  is  necesstiry  to  supplement  the  foregoing 
with  the  concept  of  an  effective  stress  (usually  the  von  Miscs  stress)  which  is  a 
function  of  a  scalar  parameter,  usually  either  the  accumulated  plastic  strain  or 
the  plastic  work.  This  approach  leads  to  the  concept  of  isotropic  hardening,  in 
which  the  yield  surface  m^ntains  its  shape  while  its  growth  in  size  depends  on 
the  pstfameter  just  cited.  Isotropic  hardening  does  not  permit  modeling  of  the 
Bauschinger  effect,  for  which  recourse  may  be  made  to  the  kinematic  hardening 
concept  of  Prager  [58]  eind  Ziegler  [59]  in  which  the  yield  surface  (in  principal 
stress  space)  retains  its  initial  size  and  shape  but  translates  a  the  direction  of  the 
plastic  deformation-increment  vector.  The  coordinates  of  the  center  of  the  yield 
surface,  usually  termed  the  back  stress  tensor,  must  be  determined  by  solution  of 
evolution  equations.  Experiment2d  evidence,  Hu  et  al.  [60],  indicates  that  hard¬ 
ening  is  neither  isotropic  nor  kinematic.  Hodge  [61]  has  discussed  a  combination 
of  kinematic  and  isotropic  hardening,  the  ratio  being  controlled  by  a  par2uneter 
varying  from  zero  to  one;  this  has  been  incorporated  in  the  D^NA  codes.  Es¬ 
sentially  all  terminal  ballistics  codes  assume  that  the  input  uniaxial  stress-st:  ain 
curve  (or  data)  is  an  odd  function  of  str^^n  (i.e.  equal  response  in  tension  and 
compression).  This  can  usually  be  made  acceptable  for  modest  strains  by  use  of 
an  appropriate  strain  measure  (generally  logarithmic  strain).  However,  for  finite 
strains  this  anti-synunetry  is  not  true  for  many  materials;  it  would  be  desirable 
to  generalize  the  constitutive  formulation  to  permit  differing  response  in  tension 
and  compression. 

Rate  Effects  Stress  evaluation  for  materials  which  exhibit  rate- dependent 
resistance  to  deformation,  i.e.  materials  which  obey  an  identifiable  ela^ito-visco- 
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plastic  constitutive  fimction,  is  readily  handled  by  conventional  codra  which  sup¬ 
ply  the  COTiponents  of  the  rate-of-defonnation  tensor  to  the  constitutive  model 
subroutine.  An  example  of  sudi  a  function  is  the  empirical  Johhson-Cook  model 
(see  «|uation  (5)),  whidi  contains  a  term  in  the  logarithm  of  the  dimensionless 
plastic  strain  rate.  In  addition  to  physical  materiad  rate  effects,  au^tificial  viscos¬ 
ity  terms  are  frequently  employed  in  numerical  simulations  to  eliminate  discon¬ 
tinuities  at  shock  fronts  aid  to  suppress  the  development  of  spurious  hourglass 
deformation  modes.  Various  constitutive  formulations  involving  strain  rates  also 
employ  the  matericJ  time  derivative  of  the  Cauchy  stress  (e.g.  the  Praudtl-Reuss 
model);  however  this  time  derivative  is  not  frame-indifferent.  Consequently,  a 
formulation  involving  the  Cauchy  stress,  its  tinie  derivative,  and  the  time  deriva¬ 
tive  of  ai  orthogonal  rotation  matrix  was  sought  which  would  provide  the  desired 
frame-indifference.  One  such  formulation  employs  the  co-rotational  stress  rate  in¬ 
troduced  by  Jaumann  [62]  which  employs  the  material  spin  tensor  (antisymmetric 
pat  of  the  velocity  gradient)  as  the  rate-of-rotation  matrix.  Although  this  stress 
rate  is  widely  used  in  current  hydrocodes  it  has  been  found  [63,  64]  to  result  in  a 
spurious  prediction  of  the  stress  in  kinematically  hadening  (hypoelastic)  materi¬ 
als  at  finite  shea  strains.  The  Jaumann  stress  .  ate  is  not  unique  in  its  capability 
to  satisfy  the  frcime-indifference  requirement.  There  is  an  extensive  literature  con¬ 
cerning  alternate  stress  rate  formulations  and  “plastic  spin”;  in  addition  to  those 
just  mentioned,  see  papers  by  Green  and  Naghdi  [65],  Asao  and  Rice  [66],  Nemat- 
Nasser  [67],  Lee  [68],  Dafalias  [69],  and  Johnson  and  Bammann  [70].  Despite  the 
proliferation  of  ideas  concerning  this  topic,  this  writer  believes  the  identification 
of  a  stress  rate  model  which  can  be  rigorously  justified  and  which  will  provide 
agreement  with  experiments  remains  to  be  accomplished. 

Energy  balance,  Thermal  Softening  Most  hydrocodes  compute  the  plas¬ 
tic  work  at  each  cycle  for  use  in  the  equation-of-state  but  few  actually  calculate 
the  temperature  field.  An  exception  is  the  EPIC-2  code  as  modified  by  John¬ 
son  [71],  which  has  an  option  to  compute  an  incremental  temperature  change  by 
dividing  the  local  plastic  work  by  the  nodal  mass  2md  the  specific  heat.  The  tem¬ 
perature  increment  is  added  to  the  temperature  from  the  previous  cycle  and  the 
updated  temperature  is  used  in  the  Johnson-Cook  [41]  flow  law  which  contains  a 
term  for  modeling  thermal  softening.  Between  each  “mechanical”  cycle  the  EPIC- 
2  code  performs  a  finite  element  heat  conduction  calculation  to  determine  a  new 
temperature  distribution  for  use  in  the  next  cycle. 

The  DYNA3D  code  has  no  provision  for  coupled  heat  transfer  calculations  but 
does  contain  the  Johnson-Cook  constitutive  model  as  Material  Model  15  (conse¬ 
quently  its  use  is  only  appropriate  for  adiabatic  applications).  Howevra,  separate 
heat  conductions  can  be  made  using  TOPAZ  [72]  to  create  files  of  temperatures  to 
be  read  by  the  DYNA  or  NIKE  codes,  thus  permitting  thermal  stress  calculations 
to  be  performed. 

The  foregoing  code  capabilities  involve  numerous  assumptions  which  may  not 
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be  generally^^alid  and  which  could  be  avoided  by  greater  complexity  of  the  formu¬ 
lation.  For  example,  it  is  assumed  that  100%  of  the.  plastic  work  is  converted  to 
heat,'  that  the  effect  of  heating  on  material  properties  is  instantaneous,  and  that 
the  specific  heat  is  a  constant.  The  writer  is  not  aware  of  any  plasticity  model  in 
any  code  which  approaches  the  generality  of  the  formulation  given  by  Green  and 
Naghdi  [65]. 

Anisotropic  Plasticity  The  writer  is  not 'aware  of  any  hydrocodes  which 
currently  provide  treatment  of  generally  imisotropic  material  in  the  plastic  re¬ 
gime.  DYNA3D  contains  a  model  for  crushable  honeycomb  (Material  Model  26) 
.which  has  been  used  with  some  success.  Initially  isotropic  materials  may  deyelop 
anisotropy  during  plastic  deformation  (texture  hardening)  but  the  computationcil 
model  does  not  characterize  the  anisotropic  parameters.  The  current  versions  of 
the  EPIC  codes  do  not  appear  to  have  any  provision  for  anisotropy.  It  should  be 
added  that  some  of  the  hydrocode  development  groups  around  the  country  are 
working  on  including  some  aspects  of  plastic  anisotropy  in  their  codes,  but  general 
release  to  the  user  community  has  not  yet  occurred. 

Failure  Modeling  There  is  a  very  extensive  literature  on  modeling  of  mate¬ 
rial  fmlure  from  both  phenomenologiccd  tind  micromechanical  viewpoints,  but  here 
mention  will  only  be  made  of  failure  models  implemented  in  currently  employed 
codes.  The  most  widely  used  failure  model  for  metals  is  the  empirical  Johnson- 
Cook  model  |73].  This  is  a  cumulative  damage  model  which  expresses  the  strain 
to  fracture  as  a  function  of  strain,  strain  rate,  temperature,  pressure,  and  equiva¬ 
lent  stress.  Once  the  fracture  criterion  has  been  satisfied  the  material  can  support 
no  shear  or  tensile  stresses  but  can  sustain  compressive  hydrostatic  pressure.  The 
implementation  of  this  and  other  phenomenological  models  in  codes  is  usually  sup¬ 
plemented  with  criteria  based  on  maximum  principal  stress  or  minimum  pressure, 
with  options  as  to  post-failure  treatment  of  stresses.  The  micromechanically-based 
Steinberg-Cochran-Guinan-Limd  materid  model  [46,  47]  also  employs  a  damage 
accumulation  parameter;  this  is  used  to  calculate  a  strength-reduction  factor  which 
is  applied  to  strength  moduli  and  the  spall  strength. 

3.2  Current  BRL  Efforts  in  Finite  Plasticity 

Joint  BRL/MTL  Program  in  Constitutive  Modeling  The  BRL  has 
provided  computational  support  for  interpretation  of  MTL  quasi-static  torsion 
tests.  This  support  included  modeling  of  the  complete  Lindholm-type  specimen 
(which  revealed  a  torsiond  buckling  problem)  and  the  development  of  a  con- 
strdned  generic  strip  model  which  greatly  ;'educes  computer  runJtimeand  memory 
requirements.  The  DYNA3D  code  has  been  adapted  to  represent  constitutive  be¬ 
havior  of  materids  such  as  copper  which  exhibit  extremely  nonlinecir  stress-strdn 
characteristics  as  well  as  the  Bauschinger  effect.  A  study  of  discrepancy  between 
uniaodd  and  torsion  test  data  has  been  completed  [74].  Further  computationd 
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support  for  interpretation  of  MTL  torsional  Hopkinson  bar  data  for  inelastic 
beating  effects  is  planned. 

BRL  Extensions  to  Constitutive  Models  in  Hydrocodes  A  new  stress 
e\^uation  subroutine  bas  been  developed  for  the  DYNA3D  code  whidi  provides: 

1.  Cboice  of  Green-Naghdi  or  Jatimamn  stress  rate; 

2.  Mixed  kinematic/isotropic  hardening; 

3.  Input  of  tabular  stress-st’-ain  data;  . 

4.  Use  of  a  ( Jj,  J3)  yield  function; 

5.  Johnson-Cook  failure  modeling  ,and  thermal  softening.  ^ 

This  formulation  will  be  exercised  by  application  to  a  variety  of  simple  applications 
and  for  interpretation  of  data  being  generated  at  MTL. 

Externally  produced  constitutive  and  failure  models  are  on  the  agenda  for 
incorporation  into  the  BRL  DYNA3D  and  evaluation:  the  Mechanical  Thresh¬ 
old  Stress  model  [75],  the  Steinberg-Cochran-Guinsm-Lund  model,  the  Zerilli- 
Armstrong  model  [42],  the  latest  version  of  Bammann’s  model  [76,  77],  and  the 
failure  models  of  J.  Johnson  [29].  The  Nemat-Nasser/Chung  plasticity  algorithm 
[78]  will  also  be  studied  for  possible  use  in  certain  applications. 

3.3  Current  Activities  Elsewhere 

Dept,  of  Energy  The  new  Eulerian  codes  MESA  and  CTH,  have  been  devel¬ 
oped  at  the  DOE  laboratories  to  take  advantage  of  new  2nd-order  acctirate  ad- 
vection  algorithms.  New  materid  models  developed  at  the  DOE  labs  include  the 
Mecbauicd  Threshold  Stress  (MTS)  model  and  the  Steinberg-Cochran-Guinan- 
Lund  model.  Both  of  these  attempt  to  relate  macroscopic  stress-strain  to  the 
underlying  microscopic  mechanical  behavior.  Much  of  the  experimental  work  in¬ 
volves  observation  of  wave  profiles  in  plane  shock  experiments.  Work  is  being 
conducted  at  the  DOE  laboratories  to  port  MESA  and  CTH  to  run  on  massively 
parallel  computers  including  the  Connection  Machine  and  the  N-Cube  Although 
massively  parallel  versions  of  these  codes  are  still  under  development,  some  of 
them  have  demonstrated  impressive  speedups  compared  to  current  vector  pro- 
erasing  machines,  when  simulating  armor/anti-armor  problems. 

URI  Center  at.  Univ.  of  California,  San  Diego  An  potentially  impor¬ 
tant  effort  at  the  USCD-URI  is  contiiined  in  the  papers  [79,  80,  78]  in  which  a 
new  algorithm  for  stress  increment  evaluation  during  plastic  loading  is  discussed. 
This  algorithm  takes  advantage  of  the  fact  that  the  plastic  strain  increment  is 
usually  the  dominant  pau-t  of  the  total  strdn  increment.  The  cleumed  ad'^ntage 
for  this  rather  complex  algorithm- is  that  it  permits  use  of  large  time  steps  in  the 


temporal  integration  of  a  simulation.  However,  for  codes  employing  a  condition¬ 
ally  convergent  explicit  integration  scheme  (as  do  most  terminal  ballistics  codes) 
the  restriction  on  step  size  necessary  for  numerical  stability  vitiates  this  advan¬ 
tage.  This  may  suggest  the  need  for  less  restrictive  timestepping  schemes  than 
are  currently  in  use. 

3.4  Conclusions 

Successful  prediption  of  stress  states  and  material  failure  in  finite  strain  plaisticity 
will  require  advances  in  many  of  the  areas  cited  above,  particularly  in  modeling  of 
stress  rates,  Bauschinger  effect,  elastic  and  plastic  volumetric  effects,  and  elasto- 
thermo- visco-plastic  coupling.  Development  of  improved  modeling  in  these  areas 
must  be  guided  by  extensive  new  experimental  techniques  and  results,  especially 
for  nonproportional  loading;  intermediate  and  high  rate  loading,  plastically  in¬ 
duced  heating  and  thermal  d'  ^adation,  and  anisotropic  plasticity.  BRL  should 
be  active  in  acquisition  ane  •  valuation  of  newly  proposed  constitutive  models  such 
as  MTS,  Steinberg-Cocl  an-Guinan-Lund,  Bammann,  and  Zerilli- Armstrong  plus 
several  micromech- aically-based  failure  models  and  should  keep  fully  abreast  of 
developments  in  ALE  codes,  other  relevant  advances  in  finite  element  /  finite 
diflference  technology  and  massively  parallel  computing. 

4  Modeling  of  Material  Response  due  to  Shear 
Band  Formation  (M*  Raftenberg,  J.  Walter) 

In  modeling  adiabatic  she2u:  band  formation,  two  distinct  aspects  of  the  problem 
are  apparent:  shesu:  band  nuchation  and  subsequent  shear  band  growth.  Ad¬ 
dressing  the  second  aspect  first,  we  may  say  for  purposes  of  discussion  that  a 
shear  band  has  nucleated  when  it  has  grown  to  a  size^,  perhaps  one  to  several 
grain  diameters,  such  that  its  subsequent  growth  is  governed  by  the  macroscopic 
(continuum)  elaistic,  thermzJ,  viscorplastic  response  of  the  material  and  the  local, 
macroscopic  loading  applied  to  the  band.  In  particular,  eXtamt  experimental  re¬ 
sults  indicate  that  once  nucleated,  adiabatic  shear  bands  are  little  affected  by  and 
often  may  completely  obliterate  any  preexisting  microstructure.  Consequently, 
most  attempts  to  model  shear  band  growth  have  been  within  a  continuum  me¬ 
chanics  framework  and  have  relied  on  phenomenological  thermo- visco-plastic  flow 
laws.  For  problems  such  as  target  plugging  and,  more  generally,  to  explore  the 
post-nucleation  dynamics  of  band  growth  and  interaction  this  approach  is  feasible 
and  has  provided  important  results. 

^Like  brittle  cracks,  adiabatic  shear  bands  are  esseniially  planar  entities  whose  in-plane  size 
„may  vary  from  a  few  to  many  thousands  of  times  the  band  thickness.  -  - - 
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ojntrast  with  shear  b^md  growth,  distributed  shear  band  nucleation  clearly 
dq>ends  <m  both  the  macro-  and  nucro-mechanical  response  of  the  material.  Es¬ 
sentially  my  locJ,  microsc Je  reduction  in  material  strength  or  viscosity,  whether 
p^xisting  or  deformation-induced,  is  a  potential  shear  band  nucleation  site. 
Moreover,  the  nucleation  mechanism  itself  need  not  be  thermally  activated  in 
order  to  pve  rise  to  an  adiabatic  shear  band.  Indeed,  once  shear  flow  localization 
stattts,  a  significjmt  contribution  from  thermal  softening  is  inevitable  provid^  only 
that  deformation  continue  at  a  sufficiently  high  rate.  Soriie  postulated  nucleation 
mechanisms  include:  griun  boimdary  sliding,  softening  due  to  microvoid  forma¬ 
tion  around  inclusions  or  second-phase  peptides',  polycrystal  textural  softening, 
anisotropic  intrapanular  plastic  response,  release  of  dislocation  pile-ups  at  grain 
bdundari^,  phase  changes.  It  is  clear  that  any  complete  model  for  distributed 
shear  band  damage  most  incorporate  some  description  of  material  inhomogeneity 
at  the  micr<^cale.  It  is  worth  noting  that  the  modd  of  Curreux  ei  al.  [15,  16,  17] 
does  so  although  in  an  indirect  manner. 

The  requirements  of  large  scade  terminal  effects  simulation  demand  that  shear 
band  modeling  be  accomplished  within  a  continuum  formulation.  This  formulation 
should  incorportate  the  effects  of  relevant  micromechamics,  yet  be  applicable  at 
the  length  scale  of  computational  discretizations*  used  in  production  armor/anti¬ 
armor  simulations.  Section  3.1  will  describe  a  recent  effort  at  BRL  towards  this 
goad.  An  adternative  approach  is  that  of  Curran  et  al.  ibid,  described  in  section 
3.2.  In  this  approach -a  “microstatisticad"  description  of  shear  band  nucleation  and 
growth  is  embedded  within  a  continuum  elastic-plastic  formulation.  Ultimately, 
models  capable  of  handling  both  isolated  amd  distributed  shear  bands  are  needed. 
For  example,  a  laurge  isolated  band  (e.g.  a  plug  boundary)  might  be  treated  as  a 
singular  zone  or  surface  (i.e.  a  slide  line)  which  would  possess  a  constitutive  struc¬ 
ture  distinct  from  that  of  the  surrounding  materiad  and  determined  by  single-band 
modeling  as  discussed  below.  Distributed  sheair  band  damage  should  be  amenable 
to  statistical  mechanical  methods  similar  to  those  being  developed  to  treat  dis¬ 
tributed  crack  damage  in  brittle  materials  by  Dienes  [81,  82].  The  goal  of  such 
analysis  is  to  obtain  continuum  formulations  which  incorporate  the  relevant  mi¬ 
cromechanics  through  an  internal  state  variable  strticture.  As  with  brittle  cracks, 
the  rasential  prerequisite  for  a  statistical  treatment  is  an  adequate  imderstanding 
of  single-  and  few-band  nucleation,  growth  and  interaction  problems. 

Returning  to  modeling  of  shear  band  growth,  we  outline  some  efforts  at  BRL 
and  elsewhere.  Much  work  to  date  has  focused  on  a  one-dimensional  “model  prob¬ 
lem”  in  which  a  slab  of  finite  thickness  and  infinite  in-plane  extent  is  subjected 
to  simple  shearing,  perhaps  .superimposed  on  uniaxial  compression.  This  is  the 
simplest  geometrical  context  in  which  adiabatic  shear  localization  can  be  studied, 

•Typical  mesh  spacings  are  0.01  to  1.0  mm,  although  for  three  dimensional  simulations  larger 
spacinp  are  usually  required. 
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but  accurately  simulates  several  important  aspects  of  the  “pressure-shear”  tests 
of  Clifton  [83]  and,  somewhat  less  accurately,  the  thin-wall  tube  torsion  test  (see 
Hartley  et  al.  [84];  Lindholm  et  al  [85]).  If  the  velocity  at  the  slab  boundaries 
is  held  constant,  the  boundaries  are  thermally  insulated  and  homogeneous  initied 
conditions  are  imposed,  then  the  system  of  field  equations  has  a  spatially  homoge¬ 
neous  solution  in  which  the  mechanical  fields  are  all  functions  of  time  alone.  Since 
for  metals  the  rate  of  strain  hardening  t3rpicaily  decreases  with  increasing  tempera¬ 
ture  and  strain,  homogeneous  deformation  will  initially  yield  monotone  increasing 
shear  stress  followed  by  monotone  decreasing  shear  stress  after  thermal  softening 
doitiinates  strain  h<irdening.  The  basic  picture  of  shear  band  formation  which 
has  emerged  from  analysis  of  the  model  problem  and  from  related  experiments  is 
as  follows:  slow  initial  vairiation  of  fields  (stress,  temperature,  strain  rate,  etc.) 
associated  with  perturbations  or  inhomogeneities,  followed  by  relatively  sudden 
localization  (shear  band  formation),  and  then  slower  post-localization  variation. 
If  the  region  containing  the  shear  band  is  also  subjected  to  sufficient  pressure, 
material  failure  may  be  considerably  delayed  or  suppressed  tdtogether. 

4.1  Modeling  of  Individual  Shear  Bands 

4.1.1  Onset  of  Localization 

Loss  of  stability  of  the  homogeneous  solution  (usually  interpreted  as  infinitesi¬ 
mal  stability)  is  a  necessary  (although  not  sixfficient)  condition  for  localization 
in  th6  model  problem.  Etirly  analyses  (see  Rogers  [14])  focused  on  identifica¬ 
tion  of  the  point  on  the  homogeneous  stress  strain  (s-7)  curve  at  which  dsfd'y 
vanishes,  typically  by  using  the  flow  law  to  expand  0  =  ds  and  then  estimating 
(or  ignoring)  various  terms  on  the  right  hand  side.  Rigorous  analyses  of  Wright 
and  Walter  [86,  87]  and  Wright  [88]  show  that  0  =  ds  is  generally  not  a  suffi¬ 
cient  condition  for  growth  of  perturbations,  much  less  for  shear  bemd  formation. 
Moreover,  even  cifter  perturbation  growth  begins  (e.g.  the  strain  field  begins  to 
develop  macroscopic  inhomogeneity),  shear  band  formation  may  not  occur  until 
much  later.  This  delay  in  band  formation  has  been  observed  experimentally  by 
Marchand  and  Duffy  [89].  Stability  analysis  of  the  model  problem  is  complicated 
by  the  fact  that  the  homogeneous  solution  is  time-vairying  and  so  the  linearized 
perturbation  equations  have  time-varying  coefficients.  Many  analyses  ignore  this 
and  draw  questionable  conclusions  based  on  constant-coefficient  techniques.  To 
date,  no  rigorous  stability  analysis  exists  for  the  model  problem  which  includes  all 
the  important  macroscopic  material  behavior:  local  conversion  of  plastic  work  to 
heat,  local  heat  conduction,  thermal  softening,  strain  and  strain  rate  hardening, 
inertia.  Moreover,  because  shear  band  formation  is  a  highly  nonlinear  process, 
linear  analyses  tend  to  provide  primarily  qualitative  results.  However,  in  several 
special  cases  important  scaling  laws  have  been  obtained  (e.g.  Wright  [88])  which 
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relate  ^ious  materia  and  problem  dependent  constants  to  critical  onset  param- 
eters  such  as:  imperfection  sensitivity,  shear  b<ind  susceptibility,  critical  strain  or 
temperature  at  she^  band  formation. 

Two  three-dimensional  extensions  of  this  sort  of  analysis  have  also  b^n 
published  by  Anand  et  at.  [90,  91]  in  which  (typically)  a  planar  “shear  band^-like 
wlocity  perturbation  is  superposed  on  a  homogeneous  shearing.  The  orientation 
of  the  shear  band  is  given  by  a  normal  vector  n  whose  value  is  obtained  as  part 
of  a  linear  perturbation  analysis.  Critical  conditions  involving  the  state  of  homo¬ 
geneous  deformation,  material  parameters  and  band  orientation  are  sought  under 
which  the  perturbation  fields  grow  rather  than  decay.  More  complete  analyses 
are  needed,  especially  those  which:  include  nonlinear  aspects  of  the  problem  (see 
Wright  [92,  93]);  pertain  to  previously  unmapped  regions  of  the  parameter  space 
(m<st  anrfysis  to  date  is  based  on  a  quasi-static  approximation  of  the  field  equa¬ 
tions);  account  for  underlying  deformations  arising  from  nonproportional  loading 
and  anisotropic  (e.g.  single  crystal)  response. 

4.1.2  Localization  and  Post-localization 

The  main  objective  in  the  model  problem  during  and  following  rapid  localization  is 
to  relate  quantities  describing  the  rapidity  and  severity  of  localization  (peak  strain 
rates  and  temperatures,  rate  of  stress  collapse,  measure(s)  of  baud  width,  etc.) 
to  material  and  other  pre-locdization  parameters.  Hoivever,  during  rapid  local¬ 
ization  the  governing  equations  are  dominated  by  extremely  nonlinear  terms  and 
no  substantive  analytical  results  exist.  It  is  probable  that  some  methods  of  com¬ 
bustion  theory  (e.g.  activation  energy  asymptotics)  may  yield  significant  insight 
but  application  of  these  techniques  is  still  in  its  infauacy.  Following  localization 
at  nominad  strain  rates  appropriate  to  ballistic  impact,  complex  wave  phenomena 
occur  (primarily  due  to  the  rapid  drop  in  stress  at  the  bamd  core)  which  depend  on 
the  inertial,  elaistic,  and  thermoviscoplastic  response  of  the  material.  Indeed,  one 
of  the  most  pressing  needs  for  improved  modeling  is  flow  law  data  at  temperatures 
from  ambient  to  meiting,  strain  rates  up  to  10^s“*  and  strains  up  to  10  or  more, 
s^  Walter  [94].  Even  numerical  simulation  Cem  be  difficult,  and  most  numeri-' 
cal  treatments  to  date  have  been  exploratory  or  have  involved  limited  parameter 
studies  on  various  material  parameters  and  the  loading  rate.  The  results  of  Batra 
and  co-workers  [95,  96,  97,  98,  99]  may  be  the  most  extensive.  It  is  clear  that 
accurate,  reliable  simulation  of  single-band  and  few-band  problems  will  require 
advanced  adapth'e  numerical  methods  as  are  now  being  developed  [100]. 

Some  numerical  results  have  been  reported  concerning  band  propagation  in 
two  dimensions;  almost  nothing  has  been  reported  in  3D.  Most  efforts  have  em¬ 
phasized  nucleation  in  that  the  band  initiates  from  a  defect,  inhomogeneity  or 
comer  singularity  and  grows  until  stopping  or  reflecting  at  the  domain  boundary. 
Such  studies  are  most  relevant  to  distributed  nucleation  when  they  attempt  to 
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model  putative  conditions  at  band  nucleation  (e.g.  single  crystal  plasticity,  st  :5s« 
contration  at  a  grain  boundary  triple  point,  etc.)  It  is  also  essential  that  accurate 
numerical  techniques  be  employed;  many  extant  studies  can  be  criticized  because 
of  artificial  dependence  on  the  computational  mesh.  Another  difficulty  is  that, 
even  with  sufficient  mesh  refinement,  the  band  structure  may  rca-'h  the  domain 
boundauy  (ending  the  computation)  well  before  it  is  fully  devjlcped.  Edgewise 
band  propagation  speeds  so  obtained  may  bear  little  relation  to  tho.ie  observed 
experimentally  for  well-formed  band  structures  such  as  plug  bcundaries  [101]. 
An  effort  is  now  underway  at  the  BRL  to  simulate  edgewise  propagation  of  plug 
boundaries. 


4.2  Modeling  Large  Scale  Terminal  Effects  dne  to  Shear 
Bands 

The  previous  section  has  reviewed  instances  in  which  deformation  localization  in  a 
viscoplastic  medium  into  an  individual  shear  band  bar.  been  successfully  moa  ded 
numerically.  The  necessaury  ingredients  were  seen  to  be  the  equations  of  motion,  a 
flow  law  that  includes  explicit  thermal  softening,  and  a  spatial  discretization  fine 
enough  to  resolve  the  deformation  gradients.  As  seen  in  Figure  8,  shear  bands  in 
RHA  exhibit  widths  as  smedl  as  6  /im,  and  finite  element  edges  in  these  analyses 
have  been  at  least  an  order  of  magnitude  smaller  than  this  dimension.  Since  the 
total  domain  size  has  been  limited  to  about  1  mm*,  it  was  feasible  to  employ  such 
small  elements. 

When  numerically  simulating  a  macroscopic  problem  in  terminal  ballistics, 
such  as  target  plate  penetration  by  a  long- rod  penetrator  or  a  shaped-charge  jet, 
■monetary  cost  and  computer  storage  considerations  prohibit  the  use  of  extremely 
small  elements.  The  goal  of  explicit  modeling  of  individual  shear  bands  must  be 
replaced  by  that  of  adeiquately  introducing  into  the  simulation  the  overall  effects  of 
shear  banding.  This  generally  involves  supplementing  the  equations  of  motion  and 
flow  law  with  other  ingredients  needed  to  model  distributed  shear  band  damage. 
This  section  will  review  two  instances  of  progress  in  this  direction. 

First  is  the  modeling  by  Raftenberg  at  BRL  of  12.7  mm-thick  RHA  plate 
perforation  by  a  copper  shaped-charge  jet  produced  by  a  warhead  fired  at  15.3 
C.D.  standoff  [102].  This  work  involve  the  insertion  into  Lagrangian  hydrocode 
EPIC-2  [103]  of  an  algorithm  to  model  shear  band  effects.  The  algorithm  consists 
of  three  parts:  (i)  a  shear  band  onset  criterion,  (ii)  a  post-onset  procedure  for  stress 
reduction,  and  (iii)  the  sli define  erosion  procedure  that  is  included  in  EPIC-2  [104], 
zmd  that  is  activated  by  a  critical  equivalent  plastic  strain,  e^rode-  three  paxts 
are  applied  at  the  element  level. 

The  onset  criterion  is  based  on  a  condition  first  proposed  by  Zener  and  Hol- 
lomon  in  1944  [13],  which  states  that  shear  banding  results  from  the  local  excess 
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of  the  rate  of  thermal  softening  over  the  rate  of  work  hardening.  The  strength 
model  in  EPIC-2  is  the  Jobnson-Cook  flow  law  [41],  which  expresses  the  flow 
stress  in  tenm  of  separable  factors  involving  work  hardening,  thermal  softening, 
and  strain  rate  hardening  or  softening  as  given  by  (5).  At  a  given  time  step,  in 
each  element  that  has  not  yet  satisfied  the  onset  criterion,  that  is  currently  un¬ 
dergoing  plastic  flow,  md  that  also  imdefwent  plastic  flow  at  the  preceding  time 
step,  the  product  of  the  first  and  last  factors  is  compared  with  its  vaine  at  the 
preceding  time  step  (the  factor  involving  eP  is  not  considered).  If  this  product  has 
diminished  since  the  previous  time  step,  the  element  is  deemed  to  have  tmdergone 
shear  band  initiation.  The  element’s  equivalent  plastic  strain  at  this  time  step  is 
denoted  efnset- 

The  stress  reduction  algorithm  used  in  Raftenberg’s  shear  bmd  modeling  pro¬ 
cedure  is  discussed  next.  Deviatoric  stresses  in  an  element  that  has  satisfied  the 
shear  band  onset  criterion  are,  at  all  time  steps  thereafter,  reduced  from  their  val¬ 
ues  as  computed  without  considering  element  damage.  The  hydrostatic  stress  is 
unaltered.  The  imposed  deviatoric  stress  reduction  is  in  proportion  to  the  amount 
by  which  the  element’s  equivalent  plastic  strain  exceeds  so  that 


where  is  the  deviatoric  stress  tensor  and  sr^  is  its  value  prior  to  reduction  due 
to  damage.  The  additional  equivalent  plastic  strain  beyond  cjnwt  corresponding 
to  which  the  element’s  deviatoric  stress  is  set  to  zero,  is  introduced  as 

an  input  parameter  to  characterize  the  material.  It  is  observed  that  elements 
that  have  satisfied  the  onset  criterion  and  that  are  then  subjected  to  this  stress 
reduction  scheme  become  susceptible  to  large  deformation. 

Slideline  erosion  is  used  to  model  removal  of  fragmented  target  and  projectile 
material.  Since  this  removal  is  activated  by  a  cutoff  on  equivalent  plcistic  strain 
in  an  element  on  the  slideline,  elements  that  have  satisfied  the  shear  band  onset 
condition  and  that  thereafter  experience  large  deformation  become  likely  candi¬ 
dates  for  future  erosion.  In  this  way  the  connection  between  shear  banding  and 
fragmentation  that  was  noted  in  Section  1.3  in  the  context  of  Figures  5  and  6  is 
introduced. 

In  [102]  variation  of  the  input  parameters,  Ac^i  and  and  a  negative 

pressure  cutoff,  (used  to  model  ductile  void  formation)  affords  good  agreement 
with  experiment  in  terms  of  final  hole  geometry  and  the  time  history  of  hole 
formation.  It  remains  to  be  seen  whether  the  parameter  values  that  produce 
agreement  with  data  fot4his  particular  case  of  standoff  and  target  plate  thickness 
will  also  produce  agreement  in  other  situations. 

Workers  at  SRI  International  have  developed  a  “microstatistiral”  approach  to 
the  problem  of  representing  shear  band  effects  in  a  macroscopic  numerical  analysis 
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[105,  106].  Their  model  has  the  ambitious  goal  of  explicity  tracking,  by  means  of 
internal  state  variables,  shear  band  properties,  including  their  local  density,  size, 
and  orientations.  The  material’s  macroscopic  behavior  is  then  related  explicitly 
to  these  properties.  The  SRI  model  has  six  components: 

1.  a  “nucleation  threshold  criterion”,  which  specifies  the  conditions  under  which 
bands  nucleate; 

2.  an  assumed  size  distribution  of  shear  bands  at  the  time  of  nucleation; 

3.  .an  equation  describing  the  rate  of  band  nucleation; 

4.  an  equation  describing  the  growth  rate  of  shear  bands, 

5.  a  model  for  the  coalescence  of  shear  bands  to  form  fragment  boundaries; 

6.  the  effect  on  the  stress  tensor  of  the  shear  band  distribution. 

A  number  of  micromechanical  assumptions  are  introduced  in  the  course  of  the 
development,  and  the  model  contains  numerous  parameters  that  characterize  a 
material.  These  parameters  have  been  evaluated  for  various  steels  and  for  DU  at 
SRI  by  means  of  a  “confined  fragmenting  cylinder”  experiment. 

The  SRI  model  has  been  implemented  into  their  working  Lagrangian  hy¬ 
drocode,  C-HEMP,  [107],  and  has  been  applied  to  at  least  two  problems  of  interest 
to  terminal  ballisticians.  One  such  problem  involves  a  fragmenting  munition  [105]. 
The  second  problem  involves  perforation  of  an  RHA  plate  by  a  rod  composed  of 
4340  steel,  where  the  geometry  and  impact  velocity  of  the  rod  were  varied  [17]. 
Impact  velocities  considered  were  in  the  range  of  500  to  1500  m/s.  At  the  lower 
.velocities  predicted  shear  b£md  damage  was  confined  to  a  thin,  cylindricail  region 
propagating  ahead  of  the  penetration  interface,  consistent  with  the  plugging  mode 
of  perforation  observed  in  their  accompanying  experiments.  At  higher  velocities 
the  plugging  mode  gave  way  to  more  widely  distributed  damage,  again  in  agree¬ 
ment  with  experiment. 

More  recently,  the  SRI  model  has  been  implemented  into  the  widely-used  La¬ 
grangian  hydrocode  DYNA2D  [108,  109].  This  code  includes  an  eroding  slideline 
algorithm,  so  that  the  SRI  model  can  now  be  applied  to  penetration  problems 
for  a  longer  duration  than  was  possible  with  C-HEMP.  The  manner  in  which  the 
SRI  workers  have  integrated  or  replaced  parts  5  and  6  of  their  procedure  with  the 
eroding  slideline  algorithm  in  DYNA2D  has  not  yet  been  reported.  A  comparison 
with  experiment  of  predicted  hole  geometry  and  shear  band  locations  obtained 
with  the  SRI  version  of  DYNA2D  has  also  not  been  reported.  The  SRI  approach, 
with  its  grounding  in  the  microstatistics  of  shear  bands,  offers  the  promise  of  a 
quantitative  understanding  of  the  role  of  shear  bands  in  determining  macroscopic 
behavior.  The  need  now  is  a  validation  of  the  approach  by  the  application  of 


this  version  of  DYNA2D  to  a  wide  range  of  ballistics  problems.  These  should 
include  penetrations  by  both  shaped-charge  jets  and  long  rods  into  a  wide  range 
of  materials. 

4.3  Conclusions 

Modeling  the  nucleation  and  growth  of  individual  shear  .bands  using  continuum 
thermo-visco-plasticity  theory  his  produced  important  insights  into  the  funda¬ 
mentals  of  this  phenomenon.  In  l-D,  contintied  analytical  work  is  needed  to 
understand  the  kinetics  of  rapid  loccdization.  Simulations  can  be  used  to  inexpen¬ 
sively  explore  the  post-locaJization  regime,  perform  parameter  studies,  and  infer 
^me  of  the  consequences  of  improved  plastic  flow  laws. 

The  very  complex  problem  of  di.'^t.ributed  shear  band  nucleation  must  be  ad¬ 
dressed  in  at  least  a  2-D  context;  ultimately,  a  3-D  treatment  will  likely  be  re¬ 
quired.  Considerable  cart-  needs  to  be  taken  to  model,  to  the  extent  possible,  the 
actual  material  microstructnres  believed  to  contribute  to  nucleation.  It  is  also  to 
be  expected  that  some  of  the  more  subtle  aspects  of  thermo- visco-plasticity  mod¬ 
els  (e.g.  textural  softening)  will  play  a  crucial  role  in  quantitative  simulation  of 
nucleation.  Analytical  techniques  which  prove  useful  in  quantifying  the  rapid  lo- ' 
calization  regime  in  1-D  may  be  applicable  in  2-D  and  3-D  as  well.  Multiple  band 
interaction  studies  are  barely  into  their  mfancy  as  yet.  Nonetheless,  a  through 
understanding  of  the  kinetics  of  band  interaction  is  essential  in  order  to  construct 
a  genuine  statistical  mechanics  of  distributed  shear  band  damage  and  to  obtain 
fully  satisfactory  continuum  models  for  distributed  damage  and  fragmentation. 

The  ability  tc  simulate  the,  effects  ot  abiabatic  shear  bands  in  macroscopic 
terminal  ballistics  problems  is  a  goal  of  great  practical  importance.'  Two  reported 
approches  have  been  discussed.  The  sanctioning  of  either  as  a  predictive  tool 
awaits  further  validation  with  experiment  q.  id  the  evaluation  of  input  parzuneters 
for  a  variety  of  relevant  materials. 

5  Brittle  Failure  and  Granular  Flow  of  Ceramics 

(M.  Scheidle?,  S.  Segletes) 

Unlike  ductile  metals,  ceramics  can  undergo  extensive  microcracking  when  sub¬ 
jected  to  dynamic  compressive  loads.  This  compressive  fracture  degrades  the 
elastic  moduli,  yield  strength  and  ten.'^ile  strength  of  the  ceramic,  resulting  in 
more  extensive  fragmentation  on  the  arrival  of  tensile  reflections  from  free  bound¬ 
aries  and  material  interfrees,  ^.nd  in  the  granular  flow  and  subsequent  bulking  of 
comminuted  material  under  shearing  loads.  The  phenomenological  observations 
in  Section  1.2  and  in  Shockey  et  al  [1!‘0]  highlight  the  influence  of  these  failure 
mechanisms  on  the  penetration  process.  These  phenomena  account,  at  least  in 
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part,  for  the  lack  of  success  in  simulating  the  penetration  of  ceramic  armor  using 
traditional  brittle  failtire  models  for  metals. 

Currently,  within  the  DOD  community  (including  DOD,  DOE,  contractors 
and  universities)  there  is  substcintial  experimental  and  modeling  activity  in  the 
area  of  dynamic  behavior  of  ceramics.  The  Ceramic  Modeling  Working  Group 
(CMWG),  hosted  by  Ed  Cort  of  ATAC,  Los  Alamos  National  Laboratory,  meets 
periodically  and  provides  a  forum  for  modelers  to  discuss  ongoing  issues  of  ceramic 
modeling.  To  date,  more  effort  has  been  placed  on  the  pre-failure  behavior  and 
the  fracturing  process  ithan  on  the  the  post-failure  response  of  reramics. 

In  the  sequel,  the  processes  of  dynamic  fracture  and  the  response  of  dcimaged 
but  intact  material  are  discussed  first.  This  is  followed  by  a  discussion  of  the 
granular  flow  of  a  comminuted  ceramic. 

5.1  Dynamic  Failure  of  Ceramics 

Theoretical  and  expv’rimental  studies  have  identified  a  variety  of  possible  mecha¬ 
nisms  for  compression-induced  microcracking;  e.g.  crack  initiation  at  grain  bound¬ 
aries,  deformation  twinning,  axid  splitting,  shear  faulting,  and  void  collapse.  The 
dominant  mechanism  may  vary  with  purity,  grain  size  and  other  microstruc- 
tural  features,'  as  well  as  with  temperature,  strain  rate  and  confining  pressure; 
cf.  (Ill,  112,  113,  114,  115,  116,  117,  118).  For  both  compressive  and  tensile 
failure,  as  strain  rate  increases  the  inertial  dependence  of  crack  instability  and 
propagation  results  in  a  higher  fracture  strength  and  a  wider  range  of  cracks 
sizes  activated  nearly  simultaneously.  Thus  damage  tends  to  be  more  evenly 
distributed  and  fragment  sizes  tend  to  decrease  with  increasing  strain  r  ite;  cf. 
(116,  118,  119,  120,  121,  122).  The  extent  of  microcrack  growth  depends  s  rongly 
on  the  duration  of  the  stress  pulse  [115],  and  there  is  some  evidence  tha;  crack 
instability  is  also  dependent  on  pulse  duration  [123].  Other  aspects  of  dvneunic 
fracture  are  not  well  understood;  e.g.  the  interaction  of  growing  microcr:  ,cks  by 
means  of  stress  waves,  and  the  influence  of  inerticil  resistance  to  crack  opening 
on  the  effective  stiffness  of  microcracked  bodies  [121].  Even  in  the  static  case, 
crack  interaction  is  an  extremely  complex  problem.  For  exaimple,  recent  v  ork  by 
Kachanov  ei  al.  [124]  indicates  that  the  effect  of  a  distribution  of  microcracis  on  a 
macrocrack  may  r^lnge  from  shielding  to  amplification  and  cannot  be  modeled  by 
-  replacing  the  microcracked  material  with  an  “effective”  elastic  material  of  r  jduced 
stiffness,  as  is  often  reported  in  the  damage  mechanics  literature. 

Much  of  the  relevant  data  on  dynamic  behavior  of  ceramics  comes  from  ]  lormal 
plate  impact  experiments.,  In  some  ceramics,  compression-induced  microciacking 
is  observed  below  the  Hugoniot  elcistic  limit  (HEL),  while  in  others  it  dees  not 
occur  until  well  above  the  HEL  [115, 125].  The  measured  HEL’s  of  many  ceramics 
are  substantially  higher  tham  the  values  calculated  fronrqurisi-static  uniaxial  strain 
tests  (even  after  the  correction  for  uniaxial  strain)  [9,  126].  HEL’s  of  2-3  times  the 
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quasi-static  ^ues  have  b^n  reported  for  alumina,  with  the  differences  attributed 
to  strain  rate  dependence  of  yield  strength  by  some  authors  and  to  pressure  de- 
p«id«»ce  of  yield  by  others  [126].  Plate  impact  experiments  [9, 10, 127]  indicate 
that  b^xjud  the  HEL  the  shear  strength  of  boron  carbide  decreases,  the  sheay 
strragth  of  alumina  and  aluminum  nitride  remauns  essentially  constant,  and  the 
shear  strength  of  silicon  carbide  and  titanium  diboride  increases  with  increasing 
shear  strain.  Because  of  the  high  confining  pressures,  retention  of  shear  strength 
and  even  increases  in  shear  strength  may  be  concurrent  with  extensive  microc¬ 
racking,  as  indicated  by  the  loss  of  spall  strength  in  AD-85  alumina,  aluminum 
nitride  imd  titanium  diboride  after  sufficiently  high  shock  loading  [10,  127], 

Constitutive  models  in  production  codes  typically  consist  of  equation  of  state 
models,  elastic-plastic  models  for  the  deviatoric  stress,  simple  failure  criteria,  and 
crude  methods  for  post-failure  strength  reduction;  cf.  the  discussion  in  Section 
2.  Even  the  more  sophisticated  models,  such  as  the  Hancock-MacKenzie  (P/Y) 
model  in  HULL  and  the  Johnson-Cook  model  in  EPIC  and  DYNA,  do  not  realisti¬ 
cally  treat  the  post-failure  strength  reduction.  In  pctfticular,  strength  reduction  is 
typically  treated  as  isotropic.  An  exception  is  the  oriented  crack  model  (Material 
Type  17)  in  DYNA  [27];  when  a  principal  tensile'  stress  excedes  a  critical  value, 
the  tensile  and  shear  stresses  are  reduced  on  the  platne  normal  to  this  principal 
direction.  However,  the  rate  of  stress  reduction  has  no  physical  basis;  the  model 
simply  “reduces  the  stresses  to  zero  over  a  small  number  of  time  steps”.  None 
of  these  models  account  for  failure  under  overall  compressive  loads.  Of  course,  it 
is  possible  to  simulate  a  restricted  class  of  experiments  with  a  phenomenological 
model  which  does  not  directly  account  for  compressive  fracture.  For  example,  uni¬ 
axial  strain  plate  impact  experiments  have  been  simulated  using  a  tensile  failure 
mode!  (for  snail),  together  with  a  plasticity  model  which  includes  the  dependence 
of  3rield  strength  and  elastic  moduli  on  such  quantities  as  temperature,  pressure, 
strain  rate,  and  a  scalar  damage  par2uneter  (which  evolves  as  a  function  of  plastic 
work);  cf.  Steinberg  [128]  amd  Furlong  tt  al.  [129].®  However,  there  is  no  evidence 
that  such  models  can  also  be  successfully  applied  to  the  dynamic  behavior  of 
ceramics  under  complex  triaxial  strain  histories. 

Of  some  relevance  to  ceramics  are  the  Mohr- Coulomb  type  models  which  have 
tradll'onally  been  applied  to  geological  materials  cind  concrete;  Material  Mod¬ 
els  16  and  25  in  DYNA  [27]  and  the  Johnson- Holmquist  model  in  EPIC  [41] 
are  <rf  this  type.  These  phenomenological  models  account  (in  a  crude  way)  for 
failure  under  overall  compressive  as  well  as  tensile  stress  states;  the  Johnson- 
Holmquist  model  dso  includes  rate  dependence  of  the  feiilure  surfaice.  To  account 
for  the  strength  of  fractured  material  under  pressure,  the  DYNA  Model  16  and  the 
Jolmson-Holmquist  model  incorporate  a  second  pressure  dependent  yield  surface 

*Steinberg*8  simulation*  indicate  that  pressure  dependence  of  the  elastic  shear  modulus  needs 
to  be  included  to  accurately  predict  the  arrival  time  of  release  waves. 


for  failed  material.  However,  shear-induced  dilatancy,  due  initially  to  the  opening 
of  microcracks  and  later  to  the  granular  flow  of  comminuted  ceramic,  is  not  ac¬ 
counted  for  in  a  physically  realistic  way  in  any  of  these  models;  cf.  the- discussion 
in  the  next  section.  •  , 

None  of  the  constitutive  models  currently  incorporated  in  production  codes 
have  all  the  features  necessary  to  predict  the  response  of  ceramic  t2irgets  under 
the  complex  loading  histories  encountered  in  a  ballistic  impact.  On  the  other  hand, 
the  problem  of  determining  which  features  need  to  be  included  in  a  cerjuxuc  failure 
model  has  yet  to  be  sorted  out.  Various  brittle  material  modeling  efforts  have  been 
put  forward  at  the  7th  CMWG  meeting  (Segletes,  [130])  by  Shockey  (SRI),  Ra- 
jendran  (UDRI),  Glenn  (LLNL),  Reaugh  (LLNL),  and  Addessio  (LANL).  Failure 
models  ba.sed  on  the  nucleation.  growth  and  coalescence  of  cracks  would  appear 
to  have  a  good  chance  of  predicting  the  time-dependent  process  of  stress  wave 
induced  fracture.  One  of  the  earliest  examples  of  this  type  (for  tensile  failure 
only)  is  the  SRI  BFRACT  model  [131].  Another  example  is  the  Statistical  Crack 
Mechanics  (SCM)  model  of  Dienes  [81*  82].  This  model  permits  crack  opening 
under  tensile  loads  and  crack  growth  under  tension  and  shear.  The  growth  of 
cracks  in  sheau'  results  in  microcrack  growth  under  compressive  loads,  so  that  the 
model  directly  accounts  (at  least  in  a  crude  way)  for  compression-induced  faulure. 
Strength  dcgredation  is  based  on  an  effective  modulus  theory,  which  naturally  ac¬ 
counts  for  both  rate  and  orientation  dependence  of  the  material  strength.  Crack 
-coalescence  is  modeled  by  a  strictly  geometric  process;  mechtinical  interactions 
are  neglected.  The  BRL  is  currently  supporting  further  development  of  the  SCM 
model  and  its  implementation  in  the  3-D  HULL  code.  Some  of  the  features  of  the 
SCM  model  are  incorporated  in  the  brittle  failure  model  of  Addessio  and  Johnson 
[132]. 

One  phenomenon  which  presents  a  difficult  modeling  problem  is  the  formation 
of  radial  cracks  under  normal  impact;  cf.  Section  1.2.  Since  this  is  initially 
an  axisymmetric  problem,  the  stress  on  every  radial  plane  is  the  same  prior  to 
the  formation  of  these  cracks,  so  that  all  radial  planes  are  equally  likely  to  fail. 
The  models  discussed  above  would  thus  predict  simultaneous  failure  on  all  radial 
planes.  However,  only  a  discrete  number  of  radial  cracks  are  observed.  An  energy- 
based  fracture  analysis  (cf.  [120,  122])  might  prove  useful  here. 

5.2  Granular  Flow 

The  granular  flow  of  brittle  materids  under  impact  is  a  topic  which  currently  con¬ 
cerns  experimenters  and  code  developers  alike  in  the  ballistics  community.  'Wlien 
considering  ballistic  penetration  into  ceramics,  it  becomes  clear  ♦hat  a  significant 
percentage  of  the  penetration  process  taJces  place  into  fractiured  ceramic.  The  di¬ 
versity  of  behavior  that  surely  must  occur  under  these  circumstances  is  daunting 
to  the  modeler:  dilatancy  effects,  which  cause  the  rubble  to  bulk  in  volume  under 


the  effects  of  shearing;  Coulomb  frictionaJ  behavior  which  inhibits  rubblized  flow 
in  a  manner  directly  related  to  the  confining  pressure  of  the  rubble;  Mid  of  course 
the  granular  flow  behavior  itself,  which  we  should  have  no  reason  to  believe  follows 
the  simple  coaxisd  flow  laws  of  metallic  plasticity.  With  this  as  the  situation,  it 
seems  a  foregone  conclusion  that  the  mechanics  of  dilatancy,  friction  and  granu¬ 
lar  flow  be  incorporated  into  the  brittle  material  models  of  today’s  hydrocodes. 
Unfortunately,  the  reality  of  the  situation  is  not  so  simply  addressed. 

Experimentally,  it  was  only  early  in  1991  that  Grady  of  SNL  [133]  devi.sed  a 
uniaxi^  plate  impact  test  to  load,  unload,  then  reload  a  ceramic  specimen,  in 
hopes  of  ascertaining  the  stress-strain  behavior  of  comminuted  ceramic. 

The  small  number  of  hydrocodes  currently  possessing  an  operational  ceramic 
model  handle  the  difficult  problem  of  granular  flow  in  a  simplistic  ad  hoc  fashion. 
For  example,  the  Johnson-Holmqiiist  model  [41]  implemented  in  the  EPIC  code 
effects  a  one-time,  instantcineous,  pressure  increase  at  the  onset  of  “fracture”  (i.e. 
when  a  scalar  damage  parameter  equals  unity),  which  is  intended  to  account 
for  the  dilatancy  of  comminuted  material  m  a  result  of  granular  flow.  Scheidler 
reports  [134]  that  this  model  is  being  modified  to  avoid  the  instantaneous  stress 
transitions  which  characterized  failure  in  the  original  model  implementation.  Note 
however,  that  this  “dilatacy”  effect  in  EPIC  is  not  explicitly  related  to  shearing 
strain  or  its  rate,  as  the  classical  definition  of  dilatancy  would  lead  us  to  believe, 
but  is  rather  tied  to  the  energy  of  compressive  volumetric  deformation.  To  further 
illustrate  the  simplistic  nature  of  cufrent  ceramic  model  development,  consider 
the  presentation  of  Henninger  (of  LANL)  to  the  7th  CMWG  [130],  in  which  he 
described  his  efforts  to  put  the  Johnson-Holmquist  model  into  the  Eulerian  MESA 
code.  The  question  Mose  as  to  how  the  model’s  scalar  damage  parameter,  which 
triggers  “instant”  failure,  could  be  a  properly  advected  quantity.  Without  specijJ 
ad  hoc  provisions  made  by  Henninger  amd  his  colleagues,  the  implementation 
would  have  permitted  failed  ceramic  material  to  “unfail",  simply  by  diffusing  its 
d.image  over  several  cells,  thus  bringing  the  damage  level  below  the  threshold 
vi.!ue  for  “failure”. 

Segletes  [135]  has  dso  presented  hydrocode  modeling  efforts  that  deal  with 
the  “bulking”  (i.e.  dilatancy)  of  ceramic  rubble  upon  failure.  Yet  here  again, 
the  evolution  of  his  bulking  terms  are  not  tied  to  shearing  rate.  The  ability  of 
this  model  to  adequately  address  granular  flow  must  thus  be  called  into  question. 
Other  recent  efforts  to  model  the  ballistic  behavior  of  ceramics  with  hydrocodes 
ignore  the  mechanics  of  granular  flow  aJtogether.  An  exception  is  the  modeling 
reported  at  the  7th  CMWG  meeting  by  Shockey  (the  model  is  actually  that  of 
Ct.“an,  Seamam  and  Cooper  [136]). 

The  granular  flow  problem  is,  by  no  meams,  an  easy  one  to  resolve.  One 
obvious  difficulty  is  the  need  to  model  the  essentiadly  discontinuous  process  of 
granular  flow  in  the  continuum  framework  6f  the  hydrocode.  Additionally,  existing 
treatments  of  low  pressure  granular  flow  (e.g.  in  the  soil  mechanics  literature). 


which  might  otherwise  provide  guidance  for  the  the  ceramic  model  developers, 
tend  to  be  characterized  by  a  large  and  “rapidly  increasing"  number  of  material 
parameters  “which  defy  physical  intuition*  (Scott,  [137]).  On  the  experimental 
side,  it  is  difficult  to  create  controlled  ^anular  flow  conditions  under  the  pressures 
and  strain  rates  characteristic  of  ballistic  impact. 

Micro-medianical  modeling  of  high  pressure  granular  flow  is  being  done  by 
Nemat-Nasser  arid  Balendran  [138]  at  the  University  of  California  San  Diego 
(UCSD),  by  Jenkins  [139]  at  Cornell  University,  and  by  Curran,  Seaman  and 
Cooper  [136]  of  SRI.  The  UCSD  model  is  based  on  the  frictional  effects  of  sliding 
granules,  wWle  the  Cornell  model,  because  of  its  ancdogy  to  the  kinetic  theory 
of  gases,  focuses  on  the  collisional  behavior  of  a  granular  mass.  What  separates 
these  models  from  others  being  developed  for  and  implemented  in  hydrocodes  is 
the  f2ict  that  the  actual  micromechanics  of  interacting  granules  are  considered  in 
the  derivation  of  the  models,  while  at  the  Stime  time  the  final  form  of  the  govern¬ 
ing  equations  are  continuum  relations  amenable  to  hydrocode  implenjentation.  It 
still  remains  to  be  seen  how  easily  the  material  parameters  characterizing  these 
models  can  be  measured. 


5.3  Conclusions 

Models  for  ceramic  materials  and  their  implementation  into  hydrocodes  is  cur¬ 
rently  in  a  great  state  of  flux,  as.  noted  by  Segletes  [130]  in  his  trip  report  from 
the  7th  Ceramic  Modeling  Working  Group  meeting: 


A  variety  of  ceramic  models  are  being  proposed  and  some  of  them 
are  being  implemented  into  hydrocodes.  What  [is  striking]  was  how 
each  computational  result  presented  seemed  different  in  fundamental 
ways;  some  caused  damage  only  under  the  penetrator  while  others 
created  damage  only  upon  tensile  reflection  at  locations  distant  from 
the  penetrator;  some  created  a  relatively  small  zone  of  damage,  while  .  ’ 
others  created  vast  zones  of  debris;  some  models  propagated  damage 
slowly  with  others  advanced  it  rapidly;  some  ceramic  bulked  while 
other  did  not;  some  rubble  was  frictional  while  other  was  not;  etc. 

As  more  experiments  with  detailed  records  become  available,  the  mettle  of 
these  new  ceramic  models  can  be  put  to  the  test.  It  is  tuo  early,  at  this  time,  to  tell 
how  useful  the  simpler  treatments  of  ceramic  modeling  in  today’s  hydrocodes  will 
be,  but  it  is  anticipated  that  modeling  which  specifically  addresses  the  mechanics 
of  compressive  frccture  and  flowing  granules  might  be  necess^lry  to  provide  the 
predictive  capability  sought  by  today’s  armor  designers. 
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6  Equation  of  State  (S.  Segletes) 

The  equation  of  state  (EOS)  provides  an  analytical  tool  by  which  the  pressure- 
volume-energy  state  of  a  material  may  be  predicted  as  the  result  of  a  thermody¬ 
namic  transition  of  the  material.  In  contrast  to  tabulated  or  graphical  state  data, 
whem  hundreds  or  thousands  of  data  points  may  be  obtained  to  diauanterize  the 
state  of  a  material  over  a  range  of  states,  the  equation  of  state,  comprised  of  semi- 
empirical  functional  forms  to  which  experimental  data  are  fit,  may  require  only  a 
few  independent  parameters  to  characterize.  Correspondingly,  relatively  few  data 
points  are  typically  acquired  to  characterize  the  EOS  parcuneters.  Often,  these 
data  are  clustered  along  a  particular  thermodynamic  path  (e.g.  the  Hugoniot,  for 
high  pressure  behavior),  while  the  behavior  off  of  this  path  is  linearized  (e.g.  as 
in  the  case  of  the  Mie-Grueneisen  EOS). 

In  small  strain  models,  an  equation  of  state  is  generally  not  employed.  Rather, 
the  hydrostatic  component  of  stress  is  simply  compiited  as  the  product  of  the 
bulk  modulus  and  the  volumetric  strain.  When  the  preslur^ volume  relationship 
is  known  to  be  non-linear,  or  when  energy  becomes  coupled  with  the  pressure,  as 
in  the  case  of  impact  events,  an  equation  of  state  must  be  used.  In  hydrocodes, 
where  event  durations  to  be  modeled  aue  typically  in  the  microsecond  timeframe, 
adiabaticity  is  a  valid  and  almost  universal  assumption.  However,  for  cases  inter¬ 
mediate  to  small  strain  and  violent  impact,  neither  constant  compressibility  nor 
adiabatic  assumptions  are  warrented,  and  heat  transfer  can  become  a  significant 
term,  which  will  affect  the  pressure  calculation.  Such  intermediate  cases  are  not 
well  handled  by  today’s  hydrocodes,  which  generally  ignore  the  effects  of  heat 
transfer  altogether.  A  heat  conduction  option  was,  in  fact,  employed  in  versions 
of  EPIC2  through  the  1980’s,  but  seems  to  have  been  dropped  when  the  EPIC2 
and  EPIC3  codes  were  combined  into  a  single  (ID,  2D,  3D)  EPIC  code  [28]. 

Problems  may  arise  in  the  application  of  equations  of  state  in  hydrocodes,  not 
only  because  of  the  sparse  data  which  was  used  to  characterize  them,  but  for  a 
variety  of  other  reasons.  Segletes  [140]  points  out  a  variety  of  problems  which  zuise 
in  the  development  and  application  of  many  of  today’s  hydrocode  EOS  models. 
They  include: 

1.  Pemutting  the  use  of  an  EOS  model  in  a  ma-lqer  which  violates  the  limita¬ 

tions  and/or  as  umptions  of  the  model.  Examples  of  this  problem  include 
EPIC  [28]  and  HUI  L  [24, 141, 142, 143],  which  do  not  limit  the  compression 
range  over  which  an  EOS  model  may  be  used  to  those  compressions  where 
the  data  fit  is  valid.  , 

2.  Relying  on  h'  -^listic  material  models  to  the  extent  where  the  models  become 
institutionalized  arid  taken  as  fact  by  an  unknowing  user  community.  Ex¬ 
amples  of  this  are  the  porous  material  models  found  in  several  hydrocodes. 
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The  HULL  code  manual  does  a  good  job  of  describing  the  ad  hoc  nature 
of  its  geological  model,  going  to  the  extent  of  warning  the  user  about  using 
this  Hugoniot  model  in  various  stress  regimes.  The  EPIC  manual  [28]  refers 
to  a  version  of  the  HULL  model  for  EPIC’s  “crushable  solid”  model,  but 
simply  describes  the  mathematics  of  the  model,  and  not  its  heuristic  founda¬ 
tions.  The  CTH  code  [144]  uses  a  porous  material  model,  and  even  assigns 
it  the  appellation  “Snowplow  model”.  No  reference  or  technical  description 
is  given,  though  a  warning  about  not  using  it  for  partially  crushed  states  is 
provided. 

3.  Making  ad  hoc  modifications  to  older,  well  founded  EOS  models,  rather 
than  rederivation  from  basic  principles.  For  example,  tk<»  EPIC  code  uses 
the  same  set  of  pressure-energy  equations  to  handle  expansion  as  well  as 
compression,  and  thus  makes  the  (presumably)  unintended  assumption  of  a 
tensile  shock. 

4.  Dissemination  of  undocumented  materi£j  libraries.  HULL  and,  to  a  lesser  ' 
extent,  MESA  are  offenders  in  this  regcird.  Such  actions  promote  the  unin¬ 
formed  use  of  data  of  unknown  origin. 

In  addition  to  these  problems,  EOS  models  which  reference  all  thermodynamic 
states  back  to  a  Hugoniot  reference  curve  (e.g.  the  Mie-Grueneisen  EOS),  can 
suffer  from  problerm  which  depend  on  the  impact  velocity  regime  of  application. 
For  low  impact  velocities,  use  of  the  shock  Hugoniot  reference  curve  is  not  correct, 
since  conditions  will  tend  to  be  isothermal,  rather  than  adiabatic.  As  contact 
velocity  increases,  plastic  work  will  increase,  tending  to  move  the  thermodynamic 
state  further  from  the  Hugoniot  reference  curve.  Since  the  Mie-Gru'eneisen  EOS 
is  linearized  about  this  reference,  the  accuracy  in  modeling  these  non-Hugoniot 
states  is  reduced.  In  the  hypervelocity  impact  regime,  additional  issues  need  to  be 
considered.  Material  phase  changes  may  play  a  significeint  role  in  the  event  but  are 
not  addressed  by  the  Mie-Grueneisen  EOS.  Fincdly,  Segletes  [140]  has  showm  that 
improper  formulation  of  the  Grueneisen  pairameter  can,  in  cases  of  hypervelocity 
impact,  overdrive  the  EOS  into  inaccurate  or  even  imstable  regimes.  All  codes 
in  his  study  (i.e.  EPIC,  HULL,  DYNA,  MESA,  and  CALE)  were  shown  to  be 
susceptible  to  this  mode  of  instability. 

The  talk  here  has  concentrated  on  the  Mie-Grueneisen  EOS,  because  of  its 
pervasive  use  in  the  hydrocode  community.  Other  options  in  equation  of  state 
modeling  axe  however,  available.  Probably  the  second  most  popular  form  in  hy- 
drocodes  is  the  Tillotson  equation  of  state  [145]  which  is  generally  applicable  at 
much  higher  pressures  than  is  Mie-Grueneisen.  Tillotson  provides  a  form  which 
essentially  allows  the  Grueneisen  coefficient  to  be  a  function  of  internal  energy, 
as  well  as  volume.  In  this  manner,  behavior  at  lower  pressures  can  be  smoothly 
transitioned  to  Thomas- Fermi-Dirac  statistical  theories  at  higher  pressures. _ 
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The  hmd^ental  undeist^diog  of  equations  of  state  are  veiy  good,  in  general. 

However,  as  appli^l  in  today’s  hydrocodes,  several  admonitions  are  applicable: 

1.  Ad  hoc  material  models  ssve  a  need,  but  are  no  substitute  for  a  properly 
founded  model. 

2.  Char^terizing  a  complicated  material  behavior  in  terms  of  a  small  handful 
of  material  parameters  can  lead  to  thermodynamic  inaccuracies  and  even 
notabilities  in  ^rious  locations  in  thermodynamic  space. 

3.  As  with  all  models,  users  must  understand  and  adhere  to  limitations  and 
assumptions  of  EOS  models  as  well  as  materisd  data.  These  limitatipns  and 
assumptions  are  not  cilways  known  to  the  code  user  community,  since  they 
are  often  ignored  or  glossed  over  in  code  documentation, 

4.  Code  developers  have  the  responsibility  to  document  material  properties,  if 
provided  with  the  code. 

5.  Correspondingly,  code  users  have  the  responsibility  to  understand  and  know 
the  origin  of  the  material  inputs  used  in  the  code. 


Summary  of  Important  Points  and 
Pressing  Needs  (T.  Wright) 


For  the  most  part,  the  comments  in  tj 
the  appropriate  place  in  the  body  of  | 
cases  may  reflect  only  the  particular  \ 
any  case,  the  intent  of  this  section  is 
of  the  mun  discussion.  These  are  noti 
point  is  followed  by  a  few  words  of  jus 
of  tasks  that  should  go  far  toward  advs 


his  section  are  amplified  and  illustrated  at 
the  report,  although  the  emphasis  in  some 
«ews  and  biases  of  the  present  writer.  In 
to  highlight  particularly  important  points 
:d  in  the  boldface  paragraphs  below.  Each 
tification  and  then  by  a  short  bulleted  list 
incing  the  state-of-the-art  in  one  particular 


^^Substantial  progress  will  proceed  from  accurate  simulation  of  carefully 
controlled  experiments.” 

This  statement,  tdcen  from  the  summary  of  Section  2.2,  was  intended  to  apply 
to  the  mechanics  of  penetration  of  cer  unic  targets,  but  it  applies  equally  well  to 
penetration  of  metallic  targets  ,^d  to  matericd  characterization  tests  of  both  ce¬ 
ramics  Mid  metaG’,  Thus  it  may  be  taken  as  a  kind  of  credo  underlying  this  whole 
report.  There  is  a  hidden  requirement,  however.  Namely,  before  it  is  pcssible  to 
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achieve  ^ctirate  numerical  simulation  of  a  complex  experiment,  it  is  first  neces¬ 
sary  to  model  the  material  response  accurately.  It  is  usually  assumed  that  the 
pointwise  response  of  a  material  is  the  same  as  its  response  to  homogeneous  de¬ 
formation.  Consequently  the  simplest  characterization  experiments  are  designed 
to  achieve  essentially  homogeneous  deformation  over  a  gauge  section  that  is  large 
enough  to  provide  access  for  measurements.  When  this  is  achieved,  there  is  no  real 
need  for  simulation.  However,  because  of  wave  propagation,  rate  effects,  complex 
test  sequences,  and  the  difficulty  of  enforcing  (pointwise)  boundary  conditions  for 
complex  specimen  geometries,  it  is  often  not  possible  to  achieve  the  degree  of  ho- 
*  ogcneity  necessary  for  direct  interpretation  of  measurements.  This  is  especially 
true  for  experiments  involving  finite  and/or  high-rate  deformations;  homogeneity 
is  definitely  not  available  as  a  simplifying  construct  in  ballistic  testing.  Under 
these  circumstances  simulation  is  essential  for  all  but  the  crudest  interpretation 
of  experiment.  •  —  -  -  -  -  - 

•  There  is  a  need  for  accurate  material  models  for  use  in  numerical  simulation 
of  ballistic  events.  In  many  cases  numerical  simulation  is  also  needed  to  help 
interpret  high  rate  material  tests. 

Material  modeling  in  ballistics  must  start  with  observed  phenome¬ 
nology. 

Tffis  statement  may  seem  obvious,  but  it  cannot  be  overemphasized.  The 
violent  physic<il  environment  of  penetrator  formation  and  bsillistic  impact  is  fzu: 
from  the  usual  training  and  experience  of  most  engineers  and  scientists.  In  fact 
the  pressures,  strains,  and  strcun  rates  routinely  encountered  in  ballistics  (aver¬ 
ages  of  the  order  of  10  GPa,  1.0,  «ind  10‘‘s~‘  to  10®s"^  respectively,  with  much 
higher  peaks,  [1])  are  difficult  to  reproduce  in  controlled  laboratory  experiments 
for  material  characterization.  By  answering  key  phenomenological  questions  such 
as  “Where  does  the  material  go?”,  “Where,  when,  and  how  does  it  deform  and 
fail?”,  or  “What  forces  exist  in  the  material  to  make  these  things  happen?”  the 
experimental  ballistician  can  provide  the  facts  that  chadlenge  those  who  develop 
material  models. 

•  There  is  an  overwhelming  need  for  careful,  diagnostic  experiments  on  ballis¬ 
tic  phenomenology.  Time  resolved  experiments  at  fixed  locations  and  time 
sequences  of  spatially  resolved  observations  are  generally  more  useful  than 
simple  post  mortem  measurements.  Microscopy  can  be  used  to  determine 
the  material  mechanisms  of  deformation  and  failure. 

There  are  significant  differences  in  the  response  of  metal  and  ceramic 
targets. 
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Agutt  this  statement  may  seem  obvious,  but  the  differences  are  often  un- 
d^^timated  in  planning  and  interpreting  ballistic  experiments.  Because  metals 
can  meperience  large  plastic  deformations,  they  axe  not  greatly  altered  bv  the 
p^sage  of  the  initial  shock  wave,  and  a  sort  of  quasi-steady  process  zone  c^ 
form  at  the  penetrator/target  interface.  This  is  particularly  evident  for  long  rod 
psietrators.*®  On  the  other  hand,  ceramics  experience  little  or  no  plastic  defor¬ 
mation  before  degradation  sets  in  so  they  are  very  sensitive  to  the  intensity  of 
w&ve  fifonts  and  wave  inter^tions.  If  process  zones  do  form  in  ceramics,  they  may 
be  lastly  different  from  those  that  form  in  metals. 

•  There  is  a  n^d  for  simple,  but  realistic,  models  of  penetration  into  ceramics. 


The  principal  numerical  methodologies  in  use  today  for  large  scale  nu¬ 
merical  simulations  of  ballistic  impacts  each  have  characteristic  draw¬ 
backs. 

Since  Eulerian  schemes  average  the  field  quantities  over  a  computational  cell, 
boundaries  between  different  materials  tend  to  become  blurred  through  a  kind 
of  numerical  diffusion  that  is  wholly  artificial.  Lagrangian  schemes  can  preserve 
material  boundaries  by  .suitably  locating  the  grid  points,  but  easily  fail  through 
grid  tangling  or  excessive  distortion,  usually  near  stagnation  points.  Other  kinds 
of  material  discontinuities,  such  as  wave  fronts,  elastic/plastic  boundaries,  or  fail¬ 
ure  boundaries,  tend  to  become  indistinct  in  either  scheme  because  they  will  not 
generally  lie  along  predetermined  grid  lines.  Eroding  interface  algorithms  in  La¬ 
grangian  codes  achieve  significantly  longer  nm  times,  but  at  the  cost  of  modeling 
the  dominamt  physical  process  with  a  numerical  artifice.  Arbitrary  Lagrangean 
Eulerian  (ALE)  type  codes,  now  coming  onstream,  would  seem  to  have  great  po¬ 
tential  for  capturing  the  best  features  of  both  major  types  while  avoiding  their 
worst  drawbacks,  but  experience  with  these  codes  is  still  limited. 

•  There  is  a  need  for  better  algorithms  to  track  material  boundaries  in  Eule- 
nss  codes,  and  to  track  other  kinds  of  discontinuities  in  both  Eulerian  and 
Lagremgean  codes. 

•  There  is  a  need  for  algorithms  to  treat  flow  near  a  stagnation  point  in  La¬ 
grangean  codes  so  as  to  achieve  unlimited  run  times  without  sacrificing  ac¬ 
curate  physical  modeling. 

*®The  consequences  of  this  process  zone  in  the  penetrator  have  been  explored  by  Wright  and 
Frauik  [146]  ,  who  noted  that  “target  resistance”  involves  characteristics  of  the  rod  and  of  the 
specific  collision  under  consideration,  as  well  as  the  flow  stress  in  the  target.  Thus  the  concept 
of  intrinsic  target  resistance,  although  often  useful,  does  not  seem  to  hold  strictly  for  metals  any 
more  than  it  does  for  ceramics. 


42 


•  There  is  a  need  to  learn  how  to  control  and  optimize  the  grid  motion  in 
ALE-type  codes. 

•  There  is  a  need  for  adaptive  gridding  algorithms  that  dynamically  move, 
refine,  or  coarsen  the  computational  grid  in  response  to  the  deformation  so 
as  to  achieve  the  highest  resolution  in  regions  with  high  gradients  in  the 
field  quantities. 

All  computing  techniques  in  current  terminal'  effects  codes  fail  where 
highly  localized  deformation  occurs. 

Physical  processes,  such  as  the  formation  of  tensile  voids,  adiabatic  sheair 
bands,  or  extensive  microcracking,  can  produce  strain  softening  amd  subsequently 
lead  to  severely  localized  deformations  that  are  well  beyond  the  resolution  of  any 
reasonable  grid  scheme.  Examples  are  spall  planes  or  plug  boundaries,  which 
should  form  in  a  calculation  as  a  natural  consequence  of  physical  processes.  Fail¬ 
ure  to  resolve  the  deformation  adequately  in  the  neighborhood  of  the  locaJization 
always  leads  to  an  overprediction  of  fiow  stress.  In  other  words,  the  finite  grid 
size  itself  becomes  a  (nonphysical)  material  property. 

•  There  is  a  need  to  learn  how  to  simulate  the  effect  of  localization  without 
resort  to  ultiafine  gridding. 

Theories  of  finite  plastic  deformation  are  still  under  development. 

.  Any  continuum  theory  of  material  behavior  must  be  based  upon  sound  kine¬ 
matics  and  must  be  properly  invariant  under,  rigid  body  rotations.  In  addition, 
there  cire  specieJ  requirements  for  a  theory  in  order  that  it  be  suitable  for  appli¬ 
cation  to  termianl  ballistics.  It  must  be  accurate  for  finite  elastic  volume  changes 
because  of  the  large  pressures  that  may  be  achieved  under  impact,  but  it  only 
.needs  to  account  for  small  elastic  distortions  since  deviatoric  stresses  are  limited 
by  the  onset  of  plastic  flow  or  fracture  processes.  It  should  also  be  accurate  for 
large,  monotoniC  plastic  strain  and  for  one  reversal  of  loading  or  one  large  path 
ch£inge  in  stress  space",'  but  it  need  not  be  accurate  for  further  cylic  deformation. 
The  assumption  of,  plastic  incompressibility  seems  to  be  adequate  for  most  pur¬ 
poses.  -Although  most  (but  not  2ill)  structural  alloys  are  relatively  insensitive  to 
strtiin  fate,  the  rates  experienced  in'bcillisti'c  events  are  so  large  that  it  is  better  to 
include  the  effect.  Thermal  effects  are  important  both  in  regions  of  large  elastic 
volume  change  and  in  regions  with  leirge  plastic  working. 

•  There  is  a  need  for  continued  theoretical  development  of  finite  deformation 
plasticity  that  is  tailored  to  the  requirements  of  impact  problems.  Correct 

"Typical  loading  and  deformation  histories,  as  experienced  along  a  material  flow  line  during 
penetration,  have  been  calculated  by  Lin  and  Batra  [147]. 


rate  equations  and  the  effect  of  large  volume  changes  seem  especially  impor¬ 
tant.  Although  it  is  desirable  that  theories  be  moti^ted  by  micromechanical 
considerations,  it  is  only  necessary  to  achieve  macrc»copic  fidelity. 

.  •  To  achieve  widespread  acceptance  there  is  a  need  to  develop  an  extensive 

library  of  physical  properties  and  elementary  material  response  functions. 
For  economic  reasons  then,  the  theory  should  be  based  only  on  standard 
physical  properties,  and  should  require  only  macroscopic  characterization  of 
materials  as  extensive  microscopy  is  generally  too  time  consuming  and  too 
expensive  for  routine  use. 

•  There  is  a  need  for  high  rate  testing  roughly  in  the  range  of  10*s”*  to  10®s“*. 
Jump  tests  from  one  rate  to  another  should  be  useful  for  discrinunating 
between  theories.  Quasi-static  tests  are  needed  only  to  establish  a  baseline. 
Testing  at  all  rates  should  extend  to  strains  of  order  1.0  or  larger.  This  can 
often  be  obtained  in  a  sequence  of  tests,  if  not  in  a  single  test. 

•  There  is  a  need  for  testing  in  «in  extended  temperature  range,  say  from  —100“ 
C  to  2/3  the  melting  temperature  with  most  testing  taking  place  within  a 
few  hundred  degrees  of  ambient.  Tests  at  elevated  temperatures  should  use 
the  highest  heating  rates  possible  so  as  to  suppress  diffusive  effects. 

•  It  would  be  useful  to  devise  high  rate  tests  with  nonproportional  load  paths. 
It  would  be  particularly  useful  if  those  tests  simulated  a  load  path  charac¬ 
teristic  of  ballistic  impact. 

Rate  effects  in  the  material  response,  whether  in  the  plastic  flow  or 
in  the  fracture  properties  of  the  target  or  the  penetrator,  imply  that 
there  will  be  scale  effects  in  ballistic  testing. 

Either  a  penetrator  or  a  target  will  be  more  effective  in  small  scede  than  in  large 
scale  if  its  material  exhibits  a  strong  strengthening  rate  effect.  This  is  because 
local  strain  rates  can  be  expected  to  scale  as  Vq/R,  the  striking  velocity  divided  by 
the  penetrator  radius.  Thus  at  the  same  impact  speed,  strain  rates  will  be  higher 
and  loading  times  shorter  in  the  smaller  scale  experiment.  As  a  consequence  the 
material  will  appear  to  be  stronger  in  smaller  scale. 

•  There  is  a  need  for  ballistic  testing  at  different  physical  scales  when  either 
of  the  materials  in%'olved  has  a  strong  rate  affect  in  its  material  response. 
When  both  penetrator  and  target  materials  exhibit  rate  effects,  the  relative 
increase  in  dynamic  strength  is  indeterminate. 

Adiabatic  shear  banding  is  a  major  mechanism  of  damage  and  material 
erosion  in  ballistic  events. 
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Plugging  in  metallic  targets,  which  has  been  known  for  a  long  time,  and  newer 
evidence  on  penetrator  erosion  indicate  the  truth  of  this  statement.  Theoreticcd 
and  experimental  study  of  the  phenomenon  continues  as  a  major  research  thrust 
today.  It  is  now  known  from  experimental  studies  [89],  confirmed  by  computa¬ 
tional  work  [96],  and  explained  by  theoretical  work  [88,,  92]  that  a  decrease  in 
the  flow  stress  does  not  necessarily  indicate  incipient  localization  although  local¬ 
ization  can  be  expected  to  occur  somewhat  later.  Rate  effects,  heat  conduction, 
and  inertia  all  have  the  effect  of  delaying  catastrophic  loss  of  strength  for  some 
time  after  thermal  softening  has  overpowered  strain  rate  and  work  hardening  in  a 
ductile  material. 

•  There  is  a  continuing  need  to  increase  understanding  of  the  fundamental 
kinetics  of  adiabatic  shear  banding. 

•  There  is  a  need  for  developing  multidimepsional  damage  models  that  use 
kinetics  as  understood  from  more  fundamental  one  dimensional  studies. 

•  Although  there  is  a  need  for  more  data  on  rate  and  temperature  effects 
in  plastic  deformation,  it  is  probably  not  necessciry  to  obtain  data  at  the 
extreme  rates  and  temperatures  that  occur  in  the  core  of  a  fully  developed 
shear  band.  This  is  because  extreme  conditions  only  occur  after  localization, 
not  before.  Furthermore,  there  is  the  prospect  of  describing  the  effect  of  the 
shear  band  on  the  surrounding  material  without  resolving  the' full  structure 
of  the  band  itself,  much  as  is  now  done  with  shock  waves. 

•  There  is  a  need  to  develop  experimental  techniques  that,  in  a  fundamental 
way,  will  allow  comparicon  of  the  susceptibility  to  shear  bands  in  different 
materials.  Of  course  adiabatic  shear  bmds  cam  be  formed  in  many  materiaJs 
by  many  different  experimental  techniques.  The  difficulty  here  is  that  shear 
band  formation  is  an  essentiaJly  inhomogeneous  phenoihenon  that  happens 
because  of  material  instability.  Without  inhomogeneity  no  localization  can 
occur,  but  in  real  materials  and  read  loading  situations  there  are  always  inho¬ 
mogeneities  of  unknown  or  uncontrolled  magnitudes.  Consequently  current 
meaisures  of  material  susceptibility,  such  as  ambient  strain  just  prior  to  lo¬ 
calization,  have  no  vailue  for  comparing  one  material  against  another. 

The  inherently  brittle  nature  of  ceramics  presents  difficulties  in  mate¬ 
rial  modeling  for  ballistic  applications  that  are  wholly  different  from 
those  in  ductile  materials. 

Ceramics  exhibit  initiation  and  extension  of  microcracking  under  compressive 
loads  [112-119].  Since  small  fractures  take  time  to  accelerate,  the  failure  process 
is  effectively  rate  dependent.  At  slow  loading  rates,  the  weakest  point  cracks  first 
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and  then  may  dominate  subsequent  failure.  At  high  loading  rates,  many  points 
can  madi  critical  levels  before  the  weakest  point  gives  way,  so  that  mmny  fr^ture 
sito  may  be  activated.  At  very  high  loading  rates  and  very  high  pressures,  as  in 
shock  waves,  massive  fracturing  and  rubblization  may  occur,  especially  in  reftected 
waws  where  mode  conversion  takes  place.  Dne  major  consequence  of  the  rate 
dependence  of  compr^sive  failure  in  ceramics  is  that  ballistic  tests  are  not  scale 
independent. 

•  Th^e  is  a  need  to  develop  high  rate  compressive  damage  and  failure  models 
for  ceramics. 

•  hi  numerical  simulation  cf  either  characterization  tests  or  ballistic  tests  in 
ceramics  there  is  a  need  to  resolve  elastic  wave  fronts  and  interactions  be- 

•  cfiuse  there  is  little  plastic  flow  before  failure. 

•  There  is  a  need  for  material  models  for  high  speed  flow  of  a  gr2uiular  matericd. 
The  important  points  here  are  the  process  of  conversion  from  intact  material 
to  fully  comminuted  material  and  the  interaction  of  presure,  volume,  and 
shearing  in  a  high  speed  flow.  In  some  cases  the  volume  expansion  that 
accompamies  comminution  can  be  substantial. 

Even  though  there  is  a  large  amount  of  accurate  data  and  information 
available  concerning  equations  of  state,  current  codes  do  not  always 
implement  EOS  models  correctly. 

•  There  is  a  need  for  better  documentation  and  care  in  applications  in  numer¬ 
ical  simulation,  rather  than  a  need  for  new  EOS  models. ' 
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Figure  1.  Hydrodynamic  IdeaUzation  of  the  Penetration  Pm«^s 


Figure  2.  Penetration  Prediction  of  Eroding-Rod  Velocitv  Models. 


Figure?.  On  logonal  Failure  Planes  and  Radial  Cracks  in  a  Recovered  Ceramic  Target  Without 
Wave  Traps. 
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inite  Residual  Penetrator  Section  for  the  flOOl  Orientation.  Photomicroi 
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